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Abstract. In this paper, basic observational properties 
of elliptical galaxies such as the integrated spectra, the 
chemical abundances and the enhancement of a-elements 
inferred from broad-band colors and line strength indices 
Mg2 and (Fe) (and their gradients), the color-magnitude 
relation, the UV fluxes, and the mass-luminosity ratios, 
are examined in the light of current theoretical interpre- 
tations, and attention is called on several points of intcr- 

■ nal contradiction. Indeed existing models for the forma- 
tion and evolution of elliptical galaxies are not able to 
simultaneously account for all of the above observational 
features. Specifically, in the context of standard star for- 
mation in the galactic-wind driven models, that are at the 
base of present-day understanding of the color-magnitude 
relation and UV fluxes, it is difficult to explain the slope 

\ of the M/Lb versus Mb relation (tilt of the Fundamental 
Plane) and enhancement of the a-elements in the bright - 

' est elliptical galaxies. We suggest that the new initial mass 
function (IMF) by Padoan et al. (1997), which depends 
on the temperature, density, and velocity dispersion of the 
medium in which stars are formed, may constitute a break- 
through in the difficulties in question. Models of elliptical 
galaxies incorporating the new IMF (varying with time 
and position inside a galaxy) are presented and discussed 
at some extent. In brief, in a hot, rarefied medium the 
new IMF is more skewed toward the high mass end than 
in a cool, dense medium, a situation which is met passing 
from low to high mass galaxies or from the center to the 
external regions of a galaxy. As a result of the changing 
IMF, the enhancement of a-elements and tilt of the Fun- 
damental Plane arc easily explained leaving unaltered the 
interpretation of the remaining properties under exami- 
nation. Finally, some implications concerning the relative 
proportions of visible stars, collapsed remnants (baryonic 
dark matter) , and gas left over by the star forming process 
are examined. 
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1. Introduction 

The conventional picture of star formation in elliptical 
galaxies is one in which the galaxies and their stellar con- 
tent formed early in the universe and have evolved quies- 
cently ever since. This view is supported by the apparent 
uniformity of ellipticals in photometric and chemical ap- 
pearance (cf. Matteucci 1997 for a recent review) and the 
existence of scaling relations, e.g. the fundamental plane 
(cf. Bender 1997). In contrast, the close scrutiny of nearby 
ellipticals makes evident a large variety of morphological 
and kinematic peculiarities and occurrence of star forma- 
tion in a recent past (Schweizer et al. 1990, Schweizer & 
Seitzer 1992, Longhetti et al. 1997a,b). All this leads to 
a different picture in which elliptical galaxies are formed 
by mergers and/or accretion of smaller units over a time 
scale comparable to the Hubble time. Furthermore, strong 
evolution in the population of early type galaxies has been 
reported by Kauffmann et al. (1996) which has been con- 
sidered to support the hierarchical galaxy formation model 
(Kauffmann et al. 1993, Baugh et al. 1996). 

Tracing back the formation mechanism of elliptical 
galaxies from the bulk of their chemo-spectro-photometric 
properties is a cumbersome affair because studies of stel- 
lar populations in integrated light reveal only luminosity 
weighted ages and metallicities, and are ultimately un- 
able to distinguish between episodic (perhaps recurrent) 
and monolithic histories of star formation and between 
star formation histories in isolation or interaction. In any 
case, there are a number of primary observational con- 
straints that should be met by any model of galaxy for- 
mation and evolution (they will be briefly summarized 
below). As nowadays most properties of elliptical galax- 
ies have been studied with sophisticated chemo-spectro- 
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photometric models in which the dynamical process of 
galaxy formation is reduced to assuming cither the closed 
box or infall scheme and a suitable law of star formation 
(Arimoto & Yoshii 1987, 1989; Bruzual & Chariot 1993; 
Bressan et al. 1994; Worthey 1994; Einsel et al. 1995; Tan- 
talo et al. 1996; Bressan et al. 1996; Gibson 1996a, b; Gib- 
son 1997; Gibson & Matteucci 1997; Tantalo et al. 1997a, 
TCBMP97). In contrast, the highly sophisticated dynami- 
cal models of galaxy formation (Davis et al. 1992; Navarro 
& Stcinmctz 1966; Hachnclt et al. 1996a,b; Katz 1992; 
Katz & Gunn 1991; Katz et al. 1996; Steinmetz & Mueller 
1995; Navarro et al. 1996; Steinmetz 1996a, b and refer- 
ences), owing to the complexity of the problem, are still 
somewhat unable to make detailed predictions about the 
chemo-spectro-photometric properties of the galaxies in 
question. Attempts to bridge the two aspects of the same 
problem are by Theis et al. (1992 and references) and most 
recently Contardo, Steinmetz & Fritze-von Alvensleben 
(1997 in preparation) and Carraro et al. (1997). How- 
ever, even in the case of the chemo-spcctro-photomctric 
approach, extant models, which ultimately stem from the 
classical galactic wind driven model of Larson (1974), of- 
ten come to solutions that are mutually discordant, when 
trying to explain one or another of the properties of el- 
liptical galaxies. In this study we would like to propose a 
coherent scenario in which a new IMF plays an important 
role and proves to be able to remove some of the points of 
contradiction. 

The plan of the paper is as follows. Section 2 shortly 
reviews the basic observational data and their current in- 
terpretation. Section 3 presents the new IMF proposed by 
Padoan et al. (1997), which is a function of the temper- 
ature, density, and velocity dispersion of the interstellar 
medium in which stars are formed. Section 4 describes our 
prescriptions for models of elliptical galaxies, in which the 
new IMF is an essential ingredient. Section 5 discusses at 
some extent the heating-cooling balance of the gas that 
is required to derive temperature, density, and velocity 
dispersion of the IMF. Section 6 presents the reference 
models and discusses some of their general properties. Sec- 
tion 7 analyzes in detail the model results concerning the 
history of star formation and chemical enrichment, the 
fractionary masses of visible and remnant stars, and the 
variation of the IMF with age and galactic mass. Section 
8 presents the broad-band colors and magnitudes and the 
mass to light ratios of the models. In particular, it presents 
our fits of the color-magnitude relation and fundamental 
plane. Finally, a summary and some concluding, critical 
remarks are given in Section 9. 

2. Learning from observations 

2.1. The G- Dwarf Analog 

As first noticed by Bressan et al. (1994, BCF94), the 
near UV spectrum (2000 A to 3500 A) of elliptical galax- 



ies shows that the relative percentage of old metal-poor 
stars is small, thus indicating that the metallicity par- 
tition function N(Z), i.e. number of stars per metallicity 
bin, cannot be the one predicted by the closed-box scheme. 
In analogy with the G-dwarf problem in the solar vicinity 
(cf. Tinsley 1980), whose solution is found in models with 
infall (Lynden-Bell 1975), Tantalo et al. (1996, TCBF96) 
show that also in the case of elliptical galaxies the in- 
fall scheme removes the above discrepancy. Similar con- 
clusion is reached by Greggio (1996), who however favors 
the prompt enrichment alternative. 

2.2. The Chemical Abundances 

Abundances in elliptical galaxies up to now have been 
deduced either from integrated colors, spectra, or linc- 
strcngth indices (Carollo et al. 1993, Carollo & Danziger 
1994a,b, Davies et al. 1993; Schombert et al. 1993) to- 
gether with their spatial gradients. In general what is mea- 
sured is a complicated result of age and metallicity, so that 
disentangling age from metallicity effects is the primary 
task of any study aimed at evaluating the metallicity of 
these galaxies (cf. the recent review by Matteucci 1997). 
Metallicity indicators such as Mg 2 and (Fe) are the most 
used to derive information on metallicity. Although pass- 
ing from the line strength indices Mg 2 and (Fe) to chemi- 
cal abundances is not a straight process (cf. Tantalo et al. 
1997b, TBC97), arguments are given to conclude that the 
mean [Mg/Fe] ratio exceeds that of the most metal-rich 
stars in the solar vicinity by about 0.2-0.3 dcx (enhance- 
ment of a-elements; O, Mg, Si, etc.. with respect to Fe), 
and the ratio [Mg / Fe] increases with the galactic mass up 
to this value (cf. Matteucci 1994, 1997). The enhancement 
of a-elements is particularly demanding as it implies that 
a unique source of nucleosynthesis has been contributing 
to chemical enrichment (type II SN from massive stars, 
the main producers of Mg). It follows from this that in 
the standard star formation scenario, i.e. with constant 
IMF, the maximum duration of the star forming activ- 
ity should be inversely proportional to the galaxy mass 
(At SF cx Mg 1 ). 

2.3. The Color- Magnitude Relation 

Elliptical galaxies obey a mean color-magnitude relation 
(CMR): colors get redder at increasing luminosity, cf. 
Bower et al. (1992) for galaxies in the Virgo and Coma 
clusters and Schweizer & Seitzer (1992) for galaxies in 
small groups and field. The cluster CMR is particularly 
narrow, whereas the field CMR is more disperse. Long 
ago Larson (1974) postulated that the present-day CMR 
is the consequence of SN-driven galactic winds. In the clas- 
sical scenario, massive galaxies eject their gaseous content 
much later and get higher mean metallicities than less 
massive ones. This implies that AtsF oc Mg, contrary to 
the trend required by the a-enhancement problem. The 
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tightness of the cluster CMR (in the U-V color) suggests 
that most galaxies are nearly coeval with ages ranging 
from 13 to 15 Gyr (Bower et al. 1992). The CMR for field 
galaxies is likely compatible with more recent episodes of 
star formation, perhaps interactions spread over several 
Gyrs (Schweizer & Seitzer 1992, Longhetti et al. 1997a,b). 
Finally, the CMR is a mass-metallicity sequence and not 
an age sequence along which bluer galaxies are signifi- 
cantly younger than red galaxies (cf. BCF94, TCBF96, 
TCBMP97, and Kodama & Arimoto 1996). 

2.4- The Age-Metallicity Dilemma 

The stellar content of a galaxy gets redder both at in- 
creasing age and metallicity thus giving rise to the well 
known age-metallicity dilemma. The Hp and [MgFe] in- 
dices are particularly suited to cast light on their separate 
effects, because Hp is a measure of the turn-off color and 
luminosity, and age in turn, whereas [MgFe] is more sen- 
sitive to the RGB color and hence metallicity. The anal- 
ysis of existing data (Gonzales 1993) leads to the follow- 
ing provisional conclusions: (i) most galaxies seem to pos- 
sess nearly identical chemical structures, i.e. high mean 
mctallicities and narrow range of mctallicities (Bressan et 
al. 1996, BCT96; Greggio 1996); (ii) galaxies do not dis- 
tribute along the locus expected from the CMR of coeval 
old objects (BCT96). In contrast, they seem to follow a 
sequence of about constant metallicity and varying age so 
that recent episodes of star formation (bursts) have been 
suggested. The major drawback with this idea is that a 
sort of synchronization is required. Looking at the differ- 
ence in Hp and [MgFe] between the nuclear region (i.e. 
Re/8 : N) and the whole galaxy (i.e. Re/2: W), and trans- 
lating AHpNw and A[MgFe]ww into AtNW and AZnw, 
i.e. in age and metallicity difference, respectively, BCT96 
pointed out that in most galaxies the nucleus is younger 
and more metal-rich than the external regions, and sug- 
gested that At NW ^ At S F oc 1/E oc 1/M G . 

2.5. The Ultraviolet Excess 

All studied elliptical galaxies have detectable UV flux 
short-ward of about 2000 A (Burstein et al. 1988) with 
large variations from galaxy to galaxy. The intensity of 
the UV emission as measured by the (1550- V) color cor- 
relates with the Mg 2 index, the central velocity dispersion 
S c , and the luminosity (mass) of the galaxy. Finally, the 
HUT observations by Ferguson et al. (1991) and Fergu- 
son & Davidsen (1993) of the UV excess in the bulge of 
M31, in which a drop-off short- ward of about 1000 A is 
detected, indicate that the temperature of the emitting 
source must be about 25,000 K. Only a small percentage 
of the 912 < A < 1200 A flux can be coming from stars 
hotter than 30,000 K and cooler than 20,000 K. 

Most likely, the UV excess owes its origin to an old 
component that gets hot enough to power the ISED of 



a galaxy in the far UV regions. Several possible candi- 
dates are envisaged (cf. Greggio & Rcnzini 1990, BCF94 
and TCBF96). The appearance of the various types of UV 
sources is governed by several important physical factors, 
each of which is affected by a certain degree of uncertainty 
still far from being fully assessed. These are the efficiency 
of mass loss during the RGB and AGB phases, the en- 
richment law AY/AZ and finally, for the specific case of 
P-AGB stars the detailed relation between the initial and 
final mass of the stars at the end of the AGB phase. Pop- 
ular sources of UV radiation are: 

(1) the classical post asymptotic giant branch (P-AGB) 
stars (see Bruzual & Chariot 1993, Chariot & Bruzual 
1991), which are always present in the stellar mix of 
a galaxy. They cannot, however, be the sole source of 
UV flux because of their high mean temperature (about 
100,000 K) and lack of sufficient fuel (cf. Greggio & Rcn- 
zini 1990). Furthermore, they hardly explain the correla- 
tion with Mg2 and S c . 

(2) Very blue HB (VB-HB) stars of extremely low metal- 
licity (Lee 1994). These stars have effective temperatures 
hotter than about 15,000 K but much cooler than those 
of the P-AGB stars. Therefore, depending on their actual 
T e ff, they can generate ISEDs in agreement with the ob- 
servational data provided that the age is let vary from 
galaxy to galaxy (Park & Lee 1997). It must be checked, 
however, whether the percentage of low metallicity stars 
is compatible with the G-Dwarf problem above. 

(3) The H-HB and AGB-manquc stars of high metallicity 
(say Z>0.07) which are expected to be present albeit in 
small percentages in the stellar content of bulges and el- 
liptical galaxies in general (Bertelli et al. 1995). Indeed, 
these stars have effective temperatures in the right inter- 
val and generate ISEDs whose intensity drops short-ward 
of about 1000 A by the amount indicated by the obser- 
vational data. With normal mass loss and AY/AZ = 2.5 
(Pagel et al. 1992), the first H-HB and AGB manque stars 
occurs at the age of about 5.6 Gyr. This age is lowered if 
AY/AZ is higher than 2.5. 

(4) Finally, the analog of the above H-HB and AGB- 
manque stars, but generated by enhancing the mass loss 
rate during the RGB phase at increasing metallicity. These 
type of stars have been named by Dorman et al. (1993, 
1995) extremely hot HB objects (E-HB). They share 
nearly the same properties of the H-HB and AGB-manque 
stars. The main difficulty with this option is the uncer- 
tainty concerning the metallicity dependence of the mass 
loss rate (cf. Carraro et al. 1996). 

The most probable channels, i.e. items (3) and (4) 
above, require a suitable metallicity partition N(Z), in 
which high metallicity bins are populated (a few percent 
of the total are fully adequate), provide a simple expla- 
nation for the dependence of the (1550-V) color on Mg 2 
(metallicity ?), velocity dispersion E c , and luminosity of 
the parent galaxy, but in turn pose strong constraints on 
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the past history of star formation and chemical enrichment 
so that the right N(Z) is generated. 

2.6. The Iron Discrepancy 

Elliptical galaxies have extended haloes of gas emitting 
in the X-ray band. Other indications of the existence of 
such haloes come from HI kinematics, planetary nebulae, 
gravitational lcnsing, and dynamical studies (cf. Carollo 
et al. 1995 and references). Measurements of the iron con- 
tent in the X-ray emitting gas (hot inter-stellar medium, 
H-ISM) from the ASCA data indicate a very low abun- 
dance, much lower (some time up a factor of 10) than 
the iron content of the stellar populations in these galax- 
ies as derived from integrated optical spectra (Arimoto 
et al. 1997). Several causes are thoroughly examined by 
Arimoto et al. (1997) going from supernova enrichment, 
chemical evolution of galaxies and clusters of galaxies, gas 
flows, heating of intra-cluster medium to reliability of the 
diagnostic used to derive iron abundances from X-ray ob- 
servations. The problem is still unsettled. 

2.7. Oxygen over- abundance and iron-luminosity relation 
of I CM 

X-ray observations of the hot intra-cluster medium (ICM) 
of galaxies by Mushotzsky (1994) show an oxygen over- 
abundance relative to iron of [O/Fe]=0.1 -j- 0.7. This result 
is an important constraint to the past chemical history of 
elliptical galaxies, since thought to be the major contribu- 
tors to the enrichment of the ICM. In brief, galactic winds 
are the favored mechanism to eject metal-rich gas (Mat- 
teucci & Vettolani 1988) provided that a natter IMF than 
the Salpeter case is used to account for the amount of Fe 
and the [O/Fe] ratio at the same time (David et al 1991; 
Mattcucci & Gibson 1995). Bimodal star formation (El- 
baz et al. 1995), while explaining the abundance ratios, 
is unable to match the CMR (Gibson 1996b). All these 
models, while recovering the total ICM iron mass, could 
only account for part of the total mass of gas measured in 
the X-ray band (Arnaud 1994). Following a suggestion by 
Trcntham (1994) that the precursors of the dwarf galaxies 
populating the faint end of the luminosity function may in 
fact be the source for 100% of the ICM gas, Gibson & Mat- 
teucci (1997) show that the bulk of the ICM gas cannot 
originate within the dwarf precursors unless the pre- wind 
binding energy was 4-5 times lower than the present day 
values. Finally, the ICM iron mass increases as a function 
of the cluster optical luminosity (E+SO galaxies) such that 
m icm _ .02iy (Arnaud ct al. 1992). 

2.8. The Fundamental Plane 

Elliptical galaxies do not populate uniformly the param- 
eter space with central velocity dispersion S c , effective 
radius R e , and surface brightness I e as coordinates. They 



cluster around a plane called the Fundamental Plane. Us- 
ing the coordinate system defined by Bender et al. (1992) 
and the identities 



L = cJeRi 
M = c 2 E 2 c R e 



(1) 
(2) 



with ci and c 2 , suitable constants, the physical coordinates 
are translated into 



kl = j=log[-] 

7 1 , r ci,M 
= ~7i ~2~L 



(3) 
(4) 
(5) 



Of particular relevance is the projection of the FP onto 
the k\ — kz plane, where the FP is seen edge on. Limited to 
the case of the Virgo elliptical galaxies to avoid distance 
uncertainties, the relation ks = 0.15fci + 0.36 with a(ks) = 
0.05 is found to hold (cf. Ciotti ct al. 1996 for details). The 
ratio (M/L) increases with galaxy mass (tilt of the FP). 

Unfortunately, current models of elliptical galaxies 
that fairly match other key properties of these systems are 
still unable to account for the tilt of the FP. Specifically, 
under standard assumptions for the IMF and star forma- 
tion rate, at any given age the models predict M/Lb ratios 
that arc nearly constant at increasing luminosity (mass) of 
the galaxy, and that increase with the age at fixed galactic 
mass. See for instance the M/Lb versus Lb relations of 
the BCT96 and TCBF96 models. An easy way out of the 
problem would be that low mass galaxies are truly young 
objects (in the sense that dominant star formation started 
much later than in less massive ones). For the arguments 
discussed above, this alternative seems to be unlikely. 

Renzini & Ciotti (1993) investigated two possible ori- 
gins of the tilt: a systematic variation of the IMF (of 
power-law type) and a trend in the relative proportions 
and distributions of bright and dark matter. The con- 
clusion is that in both cases, in order to explain the tilt 
and tightness of the FP at the same time, major changes 
and fine tuning are required. In brief, the IMF or relative 
bright/dark matter distribution should change along the 
FP, but at every position a small dispersion in the IMF 
or relative bright/dark matter distribution are required 
to preserve the tightness. In a subsequent paper along the 
same vein, Ciotti et al. (1996) looked at various effects 
of structural and dynamical nature (such as orbital ra- 
dial anisotropy, relative bright/dark matter distributions, 
shape of the light profiles) under the assumption of a con- 
stant stellar mass to light ratio. While anisotropy gives a 
marginal effect, variations in the bright/dark matter dis- 
tributions and/or shape of the light profiles can produce 
the tilt. Also in this case fine tuning is however required 
to preserve the tightness of the FP. 
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2.9. How to reconcile things ? 

It is evident that up to now no coherent explanation is 
available for the pattern of properties of elliptical galax- 
ies. Most striking is that, while the gradients in Hp and 
[MgFe] and trend for the ratio [Mg/Fe] with galaxy mass 
(Davies et al. 1993; Matteucci 1994) suggest that the du- 
ration of the star formation activity ought to be shorter 
in massive galaxies than in low mass ones, the opposite 
should occur to account for the slope of the CMR. 

In the contest of convectional picture of star formation 
in elliptical galaxies (formation and evolution in isolation), 
we address the question whether relaxing some of the stan- 
dard assumptions a viable solution can be found. As al- 
ready recalled, among the various suggestions advanced 
to explain the FP tilt, Renzini & Ciotti (1993) explored 
the effect of varying the IMF and concluded that major 
but suitable change of the IMF are required. Padoan et al. 
(1997) have proposed a new type of IMF that may turn 
out to be able to reconcile the above points of difficulty. 

3. A universal IMF 

Padoan's et al. (1997) IMF stems from a statistical de- 
scription of the density field emerging from randomly 
forced supersonic flows in star forming regions. Hydro- 
dynamical simulations of supersonic flows indicate that 
very large density contrasts can develop with log-normal 
probability distribution. The resulting IMF is 
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cj){M)dM = 



2JJ " M- 3 exp[-0.5( 2lnM - ) 2 ]dM = l. (6) 



(27TCT 2 ) 



2-10.5 ' 



The quantities A, B and a are defined by the following 
relations: 



^ = 21nB + 0.5cr 2 , 

B = 1.2 x (T/10) 1 ' 5 x (n/1000r - 5 , 

a 2 = ln[l + 0.36(A4 2 -l)], 



(7) 
(8) 
(9) 



in which a is the standard deviation of the number density 
distribution in the field with respect to the mean. M is 
the Mach number given by M 2 = (T, g /v s ) 2 , T is the tem- 
perature in K, n is the number density in cm~ 3 , v s is the 
sound velocity, and S 9 (in km/s) is the velocity disper- 
sion of the star forming gas. This IMF has a long tail at 
high masses, an exponential cutoff at the smallest masses, 
a characteristic peak mass 

M P ~ 0.2M© x (T/10) 2 x (n/T000)-°- 5 x (E g /2.5) _1 ,(10) 

and a slope continuously varying with the mass. The IMF 
gets flatter and Mp gets higher at increasing temperature 



or decreasing density and velocity dispersion. The oppo- 
site occurs reversing the variations of the three physical 
parameters. The response of the IMF to variations in tem- 
perature is shown in Fig. 1 for values of density and veloc- 
ity dispersion S g typical of molecular clouds, i.e. n — 50 
particles per cm 3 , and S 9 = 2.5 km s _1 . See Padoan et al. 
(1997) for all other details. 

The IMF is analytically integrable. Performing the 
change of variable M = e and writing 



2B 2 



a 



P = 2; 7 =A ; S= * 



(27r)°- 5 ( r' 
one gets 

/oo 
ae -/3t e -[-7(t-5) 2 ]^ 
-oo 

whose general solution is 



/* 



(/)(t)dt = 0.5[l. + Erf(x)} 



with 



x = 



f3-2j5 + 2 7 t* 



2 ( 7 )0.i 



(11) 



(12) 



(13) 



(14) 



and 

2 f x 

Erf(x) = —= I exp(-z 2 )dz. 
V 71 " Jo 

4. Modelling elliptical galaxies 

The new IMF has been implemented in a simplified ver- 
sion of the TCBMP97 models of elliptical galaxies in which 
spatial gradients of mass density (baryons and dark mat- 
ter) and star formation are taken into account. For the 
sake of simplicity and better understanding of the role 
played by the sole IMF, we adopt here the closed-box for- 
mulation (no infall of primordial gas). 

4-1. Sketch of the models 

Elliptical galaxies are modelled as a spherical distribution 
of baryonic material with total mass Ml.t (originally in 
form of gas) embedded in a large halo of dark matter with 
total mass Mp>,T- Both components have suitable radial 
distributions. For the purposes of this study, we will limit 
ourselves to consider only the region of the galaxy inside 
the effective radius R e , treated as single density- averaged 
zone. As far as dark matter is concerned this is included 
assuming that the ratios of dark to luminous mass and 
radius are M D>T /M L ^ T = ( = 5 and Rd/Rl = C = 5 (cf. 
Bertin et al. 1992; Saglia et al. 1992). Other values of ( 
are possible, whose net effect is to change the gravitational 
potential. 
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Salpeter like T=10 £ =2.5 n=50 



Fig. 1. The IMF of Padoan et al. (1997) for two different 
values of the temperature (10 and 30 K) and constant den- 
sity (n = 50 cm -1 ) and velocity dispersion E 9 = 2.5 km s _1 . 
Similar curves are possible at varying density and velocity dis- 
persion. Note how the IMF may change from a Salpeter-like 
law (very small Mp) to cases very different from the classical 
behavior. 



4-2. The basic notation 

Let Pi(r) be the density distribution of luminous matter, 
and Pi( r ) its mean value within the sphere of radius r. 
Originally this material is in form of gas and it is gradually 
turned into stars so that at each time ~p L (r, t) = p 8 (r, t) + 
~p g {r, t), with obvious meaning of the symbols. The masses 
of luminous material, gas and stars within the sphere of 
radius r are simply given by 

M L (r) =P L (r) x V(r) 

M g (r,t) =p g (r,t) x V(r) 

M s {r,t) =p s {r,t) x V(r) 

where V(r) is the volume of the sphere, supposedly con- 
stant in time. 

In the following galaxies are identified by their total lu- 
minous mass Afij, whereas the luminous mass within 
the i? e -sphere is named Ml.r^- All masses are in units of 
10 12 x M Q . 

Limited to the i? e -sphere, let us now define the a- 
dimcnsional variables 



M L,R e p L , Re 
M s (r,t) _ p s {r,t) 



G s (r,t) 



where p L R is the mean density of luminous mass within 
it. 

Finally, we introduce the gas components G 9i i(r,t) = 
G g (r 7 t) x Xi(r,t) where Xi(r,t) are the abundances by 
mass of the elemental species i. Summation of Xi(r, t) over 
all components is equal to unity. 

4-3. The spatial distribution of Ml^t o-nd Md^t 

Luminous matter. The spatial distribution of the luminous 
matter is supposed to follow the Young (1976) density pro- 
file based on the r 1 / 4 -law. Given the effective radius of a 
galaxy, one can immediately derive from the tabulations 
of Young (1976) the density pl{t), the mean density ~p L (r) 
within the sphere of radius r, and the gravitational poten- 
tial ^i(r). 

Dark matter. The mass distribution and gravitational 
potential of the dark matter are derived from the density 
profiles by Bertin et al. (1992) and Saglia et al. (1992) 
adapted to the Young formalism for the sake of inter- 
nal consistency. The details of this procedure are given 
in TCBMP97. In brief we start from the density law 



Pd(t) 



Pd,q x r^p 
{r 2 D + r 2 ) 2 



(15) 



where r^o and pd,o are two scale factors of the distribu- 
tion. 

We assume rn,a = \Rd,c, where R,D,e is the effec- 
tive radius of dark matter. Furthermore, for this lat- 
ter we assume that it scales according to the same law 
as for the total radii of dark and luminous matter, i.e. 

R]J "^Re o 6* 

The density scale factor p^o is derived from imposing 
that the total mass of dark-matter Md,t — (Ml,t, where 
Md,t is given by: 



Md,t = 47r / r 2 p{r)dr — AirpD,o r D o m (°°) (16) 
Jo 



with 

m(oo) 



T dr 



<o I 1 



\ rn,o ) 



(17) 



for which a numerical integration is required. Finally, the 
density profile of dark-matter is 



Pd(t) 



M 



D,T 



1 



1 



m(oo) 'inrjj 



1 



(18) 



The radial dependence of the dark-matter gravita- 
tional potential is 



M, 



L,R e 



Voir) = -G 



M D (r) 



dr 



(19) 
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which upon integration becomes 
VD(r) = - Air G p d, or 2 D ipb 



(20) 



where ip D (77^) is S iven b Y 
r/r D , m{;r / rDfi ) 



L 



(21) 



which must be integrated numerically. Finally, the total 
binding energy of the gas is 

n fl (r,t) = p g (r,t)V(r)ip L (r) + p g (r,t)V(r)<p D (r). (22) 

4-4- The mass-effective radius relation 

To proceed further we need to adopt a suitable relation- 
ship between the effective radius and the total mass of the 
luminous material. For the purposes of this study, we de- 
rive from the data of Carollo & Danziger (1994a, b), and 
Carollo et al. (1993, 1995) the following relation 



0.557 



R e = 17.13 x 



(23) 



with R e in kpc and M^t in units of 10 12 M Q . 
Another useful relation is between the total radius and 
total mass of the luminous material for which using the 
same sources of data we get 



(24) 



R L = 39.10 x M^j? 2 
in the same units as above. 

4-5. The chemical equations 

The chemical evolution of elemental species is governed by 
the same set of equations as in TCBF96 and TCBMP97 
however adapted to the density formalism and improved 
as far as the ejecta and the contribution from the Type la 
and Type II supernovae are concerned (cf. Portinari et al. 
1997). 

Specifically, we follow in detail the evolution of the 
abundance of thirteen chemical elements ( H, 4 He, 12 C, 
13 C, 14 iV, 16 (9, 20 Ae, 24 M ff , 28 ^, 32 S, 40 Ca, 56 Fe, and 
the isotopic neutron-rich elements nr obtained by a- 
capture on 14 A, specifically 18 0, 22 Ae, 25 Mg). Further- 
more, the stellar yields in usage here take into account 
the effects of different initial chemical compositions (cf. 
Portinari et al. 1997). 

The equations governing the time variation of the 
Gi(r,t) and hence Xi(r,t) are: 



dGi(r, t) 



dt 



-Xi(r,tMr,t)- 



/ ^(t-T M )QM, l { t ~ T M)HM)dM+ 

JM min 



,M Bm rO.5 
A 4>(M)[ f(p)*{t-TM 2 )QM,i(t-TM 2 )dp]- 

(•Bm 

{I- A) / tf (i - T M )Q M a{t - T M )<P(M)dM+ 

JMn„ 



M„ 



Mb, 



* (t - T M )Q M A t - T M )<t>(M)dM 



(25) 



where all the symbols have their usual meaning. Specifi- 
cally, ^(t) is the normalized rate of star formation, QM,i(t) 
are the restitution fractions of the elements i from stars 
of mass M (cf. Talbot & Arnctt 1973), <p(M) is the IMF, 
whose lower and upper mass limits are M min and M max 
(see below), tm is the lifetime of a star of mass M. Con- 
trary to many other models in literature (cf. Matteucci 
1997) we take into account the important variation of tm 
with the initial chemical composition using the tabulations 
by Bertelli et al. (1994). The various integrals appearing 
in eq.(25) represent the separated contributions of Type 
II and Type la supernovae as introduced by Matteucci & 
Greggio (1986). In particular, the second integral stands 
for all binary systems having the right properties to be- 
come Type la supernovae. MBm and Mbm are the lower 
and upper mass limit for the total mass of the binary 
system, f(p) is the distribution function of their mass ra- 
tios and p m i n is the minimum value of this, finally A is 
the fraction of such systems with respect to the total. We 
adopt B m = 1.5M Q , B M = 12M Q , and A = 0.2. The stel- 
lar ejecta are from Marigo et al. (1996, 1997), Portinari 
(1995), and Portinari et al. (1997) to whom we refer for 
all details. 

4.6. The star formation rate 

The rate of star formation (SFR) is assumed to depend on 
the gas density according to the Schmidt (1959) law (see 
also Larson 1991): 



dt 



= v(r,t)pJr,t) K 



(26) 



where the specific efficiency of star formation u(r,t) is a 
suitable function to be specified below. All the models we 
are going to describe are for k = 1. 
Upon normalization, the star formation rate becomes 



*(r,t) = v(r,t)[p LtR f- 1 G g (r,tr. 



4. 7. Mass limits on the IMF 



(27) 



Although the IMF of Padoan et al. (1997) is self- 
normalizing to unity over the whole range of masses from 
zero to infinity (cf. equation 6), we prefer to limit the 
mass range of validity from M min — O.O7M to M max = 



8 



IMF: new scenario 



200M Q , values that likely are more suited to real stars, 
and re-normalize the IMF by imposing that 



rM mai 
JM min 

where 
C = 



C4>{M)dM = 1 



1 - J M — cj ) (M)dM - /~ cj>(M)dM 



(28) 



(29) 



4-8. The basic energy equation 

In order to calculate the IMF we need to derive the den- 
sity, temperature and velocity dispersion S 9 of the col- 
lapsing clouds. These quantities are obtained by solving 
the energy equation 



dE(T, Pg ,t) 
alt 



■ H R (T,p g ,t) + H M + Hc 



MT, Pg ,t) 
Pa 



where E(T, p g ,t) is the energy input per unit mass, Hr, 
Hm and Hq are the heating rates (per unit mass) of ra- 
diative, mechanical, and collapse origin, respectively, and 
A is the cooling rate (in units of ergs cm -3 s _1 ). 

All physical quantities are functions of time and po- 
sition. Therefore, the IMF continuously varies with time, 
radial distance, and galaxy mass. The coupling of energy 
equation with those governing the chemical model and the 
IMF, in particular, is described in Section 5 below. 



4-9. Galactic winds 

Depending on the competition between heating and cool- 
ing, the energy stored into the interstellar medium con- 
tinuously grows thus making gas hotter and hotter. The 
case can be met in which either the thermal energy of the 
gas eventually exceeds its gravitational binding energy, or 
cquivalcntly the gas temperature exceeds the virial tem- 
perature, or the amount of gas ionized by th UV radiation 
emitted by massive stars (see below) equals the total gas 
mass. If any of these conditions is met, star formation is 
halted. The conditions in usage are 



Eth(t)>Slg{t) 



Mg tl0n > Mg 



(31) 

(32) 
(33) 



where E t h(t) and fl g (t) are the current total thermal en- 
ergy and gravitational binding energy of the gas, respec- 
tively, T g (t) is the current gas temperature, and M g ^ on 
and M g are the masses of the ionized gas and total gas, 
respectively. 



5. The energy equation: details on heating and 
cooling 

5.1. Heating 

Heating of the interstellar medium is caused by many pro- 
cesses among which we consider the thermalization of the 
energy deposit by supernova explosions (both Type I and 
II) and stellar winds from massive stars, the ultraviolet 
flux from massive stars, and finally two sources of me- 
chanical nature. 

Super novae. The rate of supernova explosions 
RsNi.n(t) over the time interval At is calculated ac- 
cording to standard prescriptions (cf. BCF94, TCBF96, 
TCBMP97, and Greggio & Renzini 1983 for Type I SN 
in particular). Knowing the amount of energy released by 
each SN explosion, the total energy injection over the time 
interval At centered at age t is the sum of the two terms 
E S ni and E S nh of type 



(30) Esniji 



j ■ 

J At 



esNi,ii(t — t')RsNi,ii{t')dt' 



(34) 



where RsNi,n{t) is the number of supernovae per unit 
galactic mass and time, and esNi,u{t) is the thermaliza- 
tion law of the energy released by a SN explosion, in which 
the effects of cooling are suitably taken into account, cf. 
BCF94, TCBF96, TCBMP97, and Gibson (1995). 

Stellar winds. The rate of energy injection by stellar 
winds is 



E 



w 



j ■ 

J At 



e w (t - t')R w (t')dt' 



(35) 



where Rw is the number of stars per unit galactic mass 
and time expelling their envelopes during the time interval 
At and, in analogy with the SN remnants, ew{t) is the 
thermalization law of the kinetic energy of stellar winds. A 
losing mass star is expected to deposit into the interstellar 
medium the energy 



(36) 



where M e j(M) is the amount of mass ejected by each star 
of mass M, v(M) is the velocity of the ejected material, 
and r\w is an efficiency factor of the order of 0.3 (Gibson 
1994 and references therein). For the velocity v(M) we 
take the maximum between the terminal velocity of the 
wind and the galaxy velocity dispersion S. For the cool- 
ing law €w(t) one could adopt the same formalism as in 
BCF94. 

Ultraviolet flux from massive stars. The rate of 
energy injection from the ultraviolet flux emitted by mas- 
sive stars, whose mass is the range 10 -j- 200 M Q , is 



E, 



uv 



J At 



e uv {t-t')R uv {t')dt' 



(37) 
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where Ruv is the number of massive stars per unit galac- 
tic mass and time and euv is the amount of ultraviolet 
energy emitted by such stars. To derive a reasonable esti- 
mate of the UV flux we adopt the following procedure. The 
problem is simplified considering that significant fluxes of 
UV radiation will be emitted only by evolutionary stages 
at high effective temperature. Therefore we limit ourselves 
to stars along the main sequence for which the relations 
between mass, luminosity, and effective temperature are 
known from stellar model calculations (cf. Bertelli ct al. 
1994). For each value of the star mass we integrate the 
spectral energy distribution up to 1000 A (our ultraviolet 
region) and calculate its ratio to the total bolometric lu- 
minosity, shortly indicated by Fuv ■ The UV flux emitted 
by each star is therefore 



(■uv = Vuv x F uv x (—)L Q 



(38) 



where (L/Lq) is the mean luminosity of a massive star 
in solar units (the mean is performed over the main se- 
quence phase), and r]uv is an efficiency parameter. The 
fraction Fuv increases from about 0.1 to 0.7 as the stellar 
mass increases from 10 to 200M Q . The effects of different 
chemical composition are neglected here. Furthermore, the 
relationship between mass and lifetime of the star is used 
to calculated the integral of eq. (37). The major prob- 
lem with the UV radiation is whether or not it must be 
included in the energy balance equation. In presence of 
dust, nearly all UV radiation is absorbed and re-emitted 
in the far infrared so that it can freely escape from the 
galaxy; indeed, the recent models of population synthe- 
sis from the far UV to the far infrared by Granato et al. 
(1997) show that about 99% of the UV radiation is re- 
processed in presence of even a modest amount of dust. 
On the basis of the above considerations, only 0.01% of 
euv above is injected into the energy equation and, ac- 
cordingly, r/uv — 0.01. Perhaps the most important use 
of the UV radiation is to set a new condition for the oc- 
currence of galactic winds. Knowing the amount of UV 
flux we calculate the radius of the Stroemgren sphere and 
associated gas mass, in which 99% of the UV flux is ab- 
sorbed and re-emitted in the infrared. When the mass of 
the Stroemgren sphere equals the total gas mass, no fur- 
ther absorption of the UV radiation is possible, so that the 
whole UV energy becomes suddenly available to the gas. 
Many numerical tests proved that, soon after, the gas gets 
so hot to exceed the virial temperature, so that condition 
(33) applies. 

Remarks on e$N and ew- According to standard 
prescriptions for esN(t), the thermalization time scale of 
SN energy is short (a few 10 6 yr) compared to the typ- 
ical time scale of chemical and photometric evolution (a 
few 10 7 to 10 8 yr). This means that only supernovae ex- 
ploding in the recent past will basically contribute to the 
current energy injection. This allows us to simplify the 



calculation of the energy input and to treat SN cooling by 
means of equation (30). The same considerations apply 
to stellar winds, for which the main source of energy are 
massive stars with short lifetimes and terminal velocities 
up to 2000 - 3000 km/s. Therefore, we drop the time de- 
pendence both in esNiji and ew- They are set equal to 
the initial energy of explosion or wind respectively. Suffi- 
ciently small time intervals At secure the desired accuracy 
and the recovery of results very close to those obtained 
with the standard treatment. 

Total radiative heating. The total heating rate 
Hji(T, p, t) of equation (30) is 



Hr(T, Pg ,t) 



Esni + Esnii + Ew + Euv 
At 



(39) 



Mechanical sources. In addition to this, we include 
two other sources of heating whose origin is of mechanical 
nature (i.e He and Hm to be defined below). 

Initial flame. The first term, He, is related to the very 
initial period of galaxy formation, when gas is collapsing 
with the infall time scale to form the galaxy. To evaluate 
this term in a simple fashion we follow Binney & Trcmainc 
(1987). Let Rh be the initial radius of the protogalaxy (say 
Rh = 100 x R e ) and r// the corresponding free-fall time 
scale 



Tff = 7T 



i? 3 
2GM 



(40) 



where M is the total mass of the galaxy, i.e. Ml.t + 
Mjj t = (1 + Q x Ml,t- In absence of energy sources of 
any type, a proto-galactic cloud at the virial temperature 
would have cooled on a time scale 



oi = 6.3 x W 5 (R h /10kpc) 2 yr. 



(41) 



Expressing Rh = XR e and using the relation between the 
effective radius and the total luminous mass, equation (23) 
above, we get 



Tcool 



120.7 x (1 + cr - 5 x A" 05 x M-° T 7S 



(42) 



with A ~ 100, and £ = 5 (for this particular set of models). 
We suppose that at the end of this phase the system has 
the energy and velocity dispersion S given by 



p — p M _ Tco ° l \ 

T ff 

E = (1 - — )°- 5 S 

T ff 



(43) 
(44) 



where E vir $ and S are the initial virial energy and ve- 
locity dispersion, respectively. The energy left over at the 
end of this phase is then converted into heat at the rate 
(per unit galactic mass and time) 



E v 



T ff 



(45) 
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This source of energy is supposed to be active as long 
as the galaxy age is smaller than r//, thus providing the 
initial gas temperature. The effect of this energy source on 
the final result is however of marginal importance, which 
somehow justifies the crudeness of the present approach. 

Cloud-cloud collisions. The second term Hm owes its 
origin to mechanical interactions of different nature among 
gas clouds, such as shock waves, friction, cosmic rays, etc. 
In order to get an estimate of this source of energy from 
simple arguments, we have looked at the cloud-cloud col- 
lisions. Limiting to the case of head-on collisions, the ki- 
netic energy per unit mass of the relative motion is £ 2 ( 
galaxy velocity dispersion), whereas the number of such 
collisions per unit time can be simply expressed by 1/tff, 
where ttt is the current free-fall time scale. Therefore the 
rate of energy injection is cast as follows 

H M = Vm x — (46) 
*// 

where t\m is an efficiency parameter masking the real com- 
plexity of the energy injection by mechanical interactions. 
Preliminary test calculations show that the parameter tjm 
ought to be of the order of 10~ 6 . This source of energy ac- 
tually plays the dominant role in the models we are going 
to describe because it will turn out to drive the minimum 
temperature attainable by the collapsing gas clouds (see 
below) . 

The CBR limit. Although not explicitly included in 
equation (30), there is another source of heating to be 
considered, i.e. the cosmic background radiation (CBR). 
It is used to limit the lowest temperature attainable by 
gas as a consequence of cooling. Specifically, we assume 
that at each time the gas temperature cannot be cooler 
than the limit Tqbr set by 

973 

T CBR [z(t)] = — x[l + z(t)} (47) 

where z(t) is the red-shift corresponding to the galactic 
age t. The red-shift z(t) depends on the choice of the cos- 
mological model of the Universe, i.e. Hq, qo (Hubble con- 
stant and deceleration parameter, respectively, in a Fried- 
mann model) , and red-shift of galaxy formation Zf or . All 
of the results below are for the standard case of Hq = 50 
km s" 1 Mpc -1 , qo = 0, and Zf or = 5. 

5.2. Cooling 

The cooling term A(T,p g ,t) of equation (30) is derived 
from literature data, and it includes several radiative 
processes. For temperatures greater than 10 4 K we lean 
on the Sutherland & Dopita (1993) tabulations for a 
plasma under equilibrium conditions and metal abun- 
dances [Fe/H]— -10 (no metals), -3, -2, -1.5, -1, -0.5, (so- 
lar), and 0.5. For temperatures in the range 100 < T < 10 4 
the dominant source of cooling is the H2 molecule getting 




Fig. 2. The net cooling rate A x n^n^ in ergs cm 3 s 1 as a 
function of temperature: from Sutherland & Dopita (1993) for 
T > 10 4 , Hollenbach & McKee (1979) and Tegmark et al. 
(1996) for 100 < T < 10 4 , and Caimmi & Secco (1986) and 
Theis et al. (1992) for T < 10 2 K. The product nj x n k con- 
tains the number densities of particles. The meaning of n-j and 
?ifc changes according to the portion of the curve under con- 
sideration: in Sutherland & Dopita (1993) nj x n k — n e X 7ij, 
where n e and n» are the number densities of electrons and sum 
of the ion number densities, respectively. In all the others, the 
product ?ij x rife = n 2 , with n the number density of all par- 
ticles lumped together for the sake of simplicity. The number 
density in usage is n ~ 10. From bottom to top, each curve 
corresponds to an increasing value of the metallicity [Fe/H]. 
See the text for more details. 

rotationally or vibrationally excited through a collision 
with an H atom or another H2 molecule and decaying 
through radiative emission. The data in use have been 
derived from the analytical expressions of Hollenbach & 
McKee (1979) and Tegmark et al. (1996). Finally, for tem- 
peratures lower than 100 K, starting from the relation of 
Theis et al. (1992) and Caimmi & Secco (1986), we have 
implemented the results of Hollenbach & McKee (1979), 
and Hollenbach (1988) for CO as the dominant coolant. 
The following analytical relation in which the mean frac- 
tionary abundance of CO is given as a function of [Fe/H], 
is found to fairly represent the normalized cooling rate (i.e. 
Aco /n- 2 with n the number density of particles) 

^ = l. 6x 10-29l ([^/ff]-l-699) T 0.5 erg cm -3 g -l (48) 

The net cooling rate (A x n 2 ) over the whole tem- 
perature range is shown in Fig. 2, in which the effects of 
metallicity are clearly visible. The number density in usage 
for the purposes of Fig. 2 is n ~ 10. It is worth mention- 
ing that no re-scaling of the cooling rate from the various 
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sources has been applied to get the smooth curves shown 
in Fig. 2. 




Gibson 1994 
Z = 0.0001 



Fig. 3. Cooling law of supernova remnants esjv(i) from Gib- 
son (1994). Each curve is for a different metallicity of the 
interstellar medium as indicated. Superposed is the cooling 
law of eq.(30) for the 3 Ml,t model with age of 0.024 Gyr 
and metallicity Z — 0.010 (thick line) described in the text. 
Re-normalization of the input energies is applied. 



How does equation (30) compare with respect to the 
classical cooling laws for supernova ejecta ? To answer 
this question we display in Fig. 3 the path in the energy- 
cooling time plane followed by supernova ejecta according 
to the most recent discussion of the subject by Gibson 
(1995) and one of the models we are going to present be- 
low, i.e. the 3 Ml,t galaxy during its stages of stellar 
activity. The cooling laws of supernova ejecta displayed 
in Fig. 3 are for different metallicities as indicated, typical 
density of the medium of about 10 -22 g/cm 3 , and injected 
energy of 10 51 erg. In the model galaxy the metallicity is 
about Z = 0.010, the parameter of H M is r\ M = 2 x 10~ 6 , 
the density is about 10 -23 g/cm 3 , and the energy injected 
is about 10 11 erg per unit mass of gas. Re-normalization 
of the energies is applied. Although the general trend is 
similar to that of supernova cooling laws, at least as far as 
the drop-off is concerned, the cooling time is much shorter. 

5.3. Coupling energy equation and IMF 

The energy equation (30) regulates heating and cooling of 
gas according to the following scheme: 

— At each time step At, the energy input from the var- 
ious sources heats up and expands the gas, whereas 
cooling acts in the opposite sense. 



Equation (30) is integrated in time assuming that 
heating-expansion and cooling-contraction of gas fol- 
low the constant-pressure law (P — -^-pT = const). 
An equilibrium stage with temperature T eq and density 
p g ,eq is soon reached that yields the Jeans mass Mj^ eq 
above which elemental gas clouds are gravitationally 
unstable and prone to collapse. 

In principle we can suppose that of the existing clouds 
only the fraction whose Jeans mass is larger than Mj^ eq 
can form stars. In order to determine the tempera- 
ture, density and velocity dispersion at the start of 
the star forming process we need to follow the thermo- 
dynamical fate of the gas component prone to instabil- 
ity. This is conceived as thermally decoupled from the 
surrounding medium, so that no further energy input 
of radiative nature is possible, but still able to mechan- 
ically interact with it. Therefore, the radiative term 
Hr(T, p g ,t) of equation (30) is switched off, whereas 
the mechanical term Hm is retained. 
The fraction of gas with Jeans mass in excess of Mj^ eq 
is let cool down and collapse until a new equilibrium 
stage is reached characterized by certain values of tem- 
perature and density, Tknee and Pknee, respectively. 
The existence of the knee-stage is secured by the en- 
ergy input of mechanical nature. In any case, Tknee is 
not let be lower than TcBR[z{t)\. In such a case the in- 
tegration of the energy equation is stopped, and Tknee 
and pknee are the last computed values. 
During all these stages the cooling and free-fall times 
scales turn out to be nearly identical. The temporal 
evolution of the gas energy (per unit mass of gas) and 
the temperatures T eq and Tknee are shown in Fig. 4 
for the SMl t galaxy at different evolutionary stages 
as indicated. The efficiency of star formation charac- 
terizing this model will be defined below. 
We identify the fraction of gas that underwent the 
cooling process described above and was able to reach 
Tknee and pknee as the molecular component in which 
fragmentation according to the Padoan et al. (1997) 
scheme can occur. 

It is reasonable to assume that these clouds acquire 
their kinetic energy from the surrounding medium un- 
der the condition of energy equipartition. Accordingly 
their velocity dispersion scales as 



E 9 (T) x 



' Mj(T) 
Mj(T eq ) 



(49) 



with obvious meaning of the symbols. The starting 
value in the recursive process to determine the veloc- 
ity dispersion of collapsing clouds is the galaxy velocity 
dispersion E. 

Finally, the temperature and density Tknee and pknee-, 
and the velocity dispersion H g (T eq ) are the values to 
be plugged into the IMF. Each cloud prone to star 
formation will fragment into stars according to this 
IMF. 
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— The case may arise, in which for the values of Tfcnee, 
Pknee, and velocity dispersion S g (T eg ) derived from 
the energy equation the peak mass of the IMF is 
comparable or even lower than the minimum mass 
M m i n = 0.07-Mq for the existence of real stars (i.e. 
burning hydrogen in their cores). In such a case the 
peak mass is not let fall below M; = O.1M , and the 
IMF is calculated with the temperature obtained from 
inverting equation (10) in which pknee and S(T eg ) are 
used. Even if there is a marginal inconsistency, this can 
be neglected for any practical purpose. 




Fig. 4. The cooling law the 3 Ml,t model with v = 10. Each 
curve corresponds to a different age and metallicity. The galac- 
tic age (in Gyr) is annotated along each curve. The particular 
stages at which the equilibrium temperature (T eq ) is reached 
and rapid contraction (collapse) begins, and the final temper- 
ature (Tknee) at which fragmentation is supposed to occur are 
indicated. See the text for more details. 



5.4- Deriving the specific efficiency of star formation 

The discussion above has been referred to the fraction 
of gas with Jeans mass in excess of Mj, eq and ready to 
collapse and fragment into stars. This fraction is, however, 
still unknown. Indeed to be determined a law governing 
the mass or number distribution of such clouds is needed. 
Similarity arguments would suggest that the same IMF 
governing star formation should also regulate the cloud 
mass (number) distribution. The major difference is in 
the temperature, density, and velocity dispersion to be 
adopted. It goes without saying that T eq and p eq must 
be used as they represent the stage at which elemental 
clouds become thermo-dynamically distinct units out of 



the surrounding medium. In such a case, the mass fraction 
of unstable clouds would be 

/•OO 

f jeans = / C 4>(M ,T eq , p g . eq , T, g (T eq ))dM (50) 

In principle, with aid of this quantity one could de- 
rive v(r, t). Recalling that cooling occurs on a certain time 
scale t coo i, known from solving equation (30), we expect 
that the rate of star formation will be inversely propor- 
tional to t coo ;. On the base of the above arguments, we 
could write the specific efficiency v(r, f) of the star forma- 
tion rate as 

/ t \ J jeans / r- -i \ 

v(r,t) = v x^- (51) 

7~cool 

where vq is a suitable scale factor of the order of unity. 

With this formulation for v{r, t), there would no longer 
be free parameters in the rate of star formation, which 
would be tightly coupled to the energy storage of the 
medium in which stars are born. Star formation drives 
the energy input, and this latter in turn affects star for- 
mation: a self-regulating mechanism is soon established. 

The application of the above scheme rises, however, 
some difficulties both of conceptual and numerical nature. 
First it is not granted that the number (mass) distribution 
of unstable clouds is governed by the same law determin- 
ing the mass distribution of stars at each formation event, 
second both f jeans and t coo i fluctuate from time step to 
time step, thus rendering the associated value of v(r, t) 
highly uncertain. The preliminary analysis of the problem 
has shown that the quantity f jeans l T coii is of the order 
of 5 4- 50. Therefore, we considered the specific efficiency 
v{r,f) as an adjustable parameter to be eventually fixed 
by comparing model results with observational data (e.g 
CMR, chemical abundances, etc..) and did not make use 
of the above scheme. 

6. The reference models 

In this section we present our reference models for ellip- 
tical galaxies of different total luminous mass Ml^t and 
luminous mass Ml^ £ within the effective radius. 

All the models are calculated assuming ijm = 2 x 10~ 6 
(the parameter governing the efficiency of mechanical en- 
ergy sources) and v — 10 (the parameter driving the effi- 
ciency of the star formation rate), but two slightly differ- 
ent conditions for the onset of the galactic wind. The main 
reason for it resides in the fact that the energy equation is 
solved separately from those governing gas consumption, 
star formation, and chemical enrichment. Therefore, there 
is an implicit ambiguity concerning the mode and time at 
which the conditions for galactic winds are met. In short, 
what is uncertain is whether the energy liberated by the 
various sources going into the thermal content of the gas 
ought to be tested against wind conditions before of after 
cooling has taken place. In order to check the sensitivity 
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of model results to this kind of uncertainty, at each time 
step we examine two possible schemes: 

— At the beginning of the time step, all the energy gen- 
erated by supernova explosions, stellar winds, etc. is 
tested against the condition 

E th {t) > n g (t). (52) 

In the case of a negative answer (no galactic wind yet), 
the total energy budget is given to the energy equation 
and all remaining conditions for the onset of galactic 
wind are examined only at the end of the time step. 
These models are named Type A. 

— Alternatively, at the beginning of the time step the 
total energy budget is given to the energy equation 
and the whole set of conditions for the onset of galactic 
winds is tested only at the end of the time step. These 
models are named Type B. 

The two sets of models are presented in Table 1, which 
lists the total luminous mass Ml^t, the total luminous 
mass Ml,r b inside the effective radius, the mean mass 
density within R e , the galaxy velocity dispersion E, the 
maximum value for the gas temperature T g (in K) and 
peak mass Mp (in M©) of the IMF reached in the course 
of galactic evolution, the age t w (in Gyr) at the onset of the 
galactic wind, the fractionary masses in gas, M 9 /Ml,r s , 
and stars M s /Mt,r c inside the i? e -sphere at the onset 
of the wind, the current (Z) and mean metallicity {(Z)), 
and finally the index Mg^ at the same epoch. The whole 
discussion below will refer to models with tjm = 2 x 10~ 6 . 
However, when examining the CMR and FP (section 8 
below), we will also consider models of Type A but with 
Vm = 3x 10 -6 which are much similar to those described 
here. 

6.1. General properties 

All the models obey some general properties that are de- 
termined by the interplay between the physical conditions 
of the star forming gas and the IMF. We will sketch the 
problem in several steps for the sake of an easier under- 
standing: 

— At increasing galaxy mass the mean gas density gets 
lower, the velocity dispersion gets higher together with 
the initial gas temperature (the one determined by the 
very first energy source He). 

— As a consequence of it, the IMF tends to skew towards 
the high mass end (Mp moves to higher masses and the 
slope above it gets slightly flatter at increasing galaxy 
mass). This means that more energy from supernova 
explosions and stellar winds is injected into the inter- 
stellar gas making it hotter and hotter till a sort of 
equilibrium stage is reached, during which both the 
gas temperature and Mp remain nearly constant. 

— However, at proceeding star formation, the gas gets 
richer in metals thus increasing the cooling efficiency 
A situation is eventually met in which cooling is so 



3 u=10 




Fig. 5. The gas temperature (in K) as a function of the age (in 
Gyr) for the models galaxies with Ml,t equal to 0.001, 0.01, 
0.1, 1 and 3 as indicated. The same curves apply both to Type 
A and B models. 




M ml„ 



Fig. 6. The peak mass Mp (in solar units) of the IMF as a 
function of the age (in Gyr) for the models galaxies with Ml,t 
equal to 0.001, 0.01, 0.1, 1 and 3 as indicated. The same curves 
apply to both Type A and B models. 
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Table 1. Basic data of the models. Ml,t and Ml,r c are the total luminous mass and the mass within the effective radius, 
respectively. They are in units of 10 12 x A/q. (p) is the mean mass density (in gr cm -3 ) within the i? e -sphere. £ is the galaxy 
velocity dispersion. T g and Mp are the maximum values of the gas temperature (in K) and peak mass (in Mq) of the IMF in 
the course of galactic evolution. t w is the age (in Gyr) of the galaxy at the onset of galactic wind. M 3 /Ml,r s and Ms /Ml, r b are 
the fractionary masses of gas and stars, respectively, at the galactic wind stage. Z and (Z) are the current and mean metallicity 
reached by the galaxy at the wind stage. All the models are calculated with rjM = 2 x 1CP 6 . Finally, two groups of models 
are listed: models of case A have early galactic wind; models of case B eject their wind at ages that progressively increase at 
decreasing galactic mass. Those with t w = 16.45 continue to form stars albeit at minimal levels up to the present epoch. 
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Threshold Metallicities 



efficient that the gas temperature and peak mass of 
the IMF in turn start decreasing. This implies that less 
energy is injected into the interstellar medium with 
consequent catastrophic fall-off of both temperature 
and peak mass. 

The temporal evolution of the gas temperature, and 
peak mass is shown in Fig. 5 and Fig. 6, respectively, 
whereas the relationship between Mp and metallicity is 
presented in Fig. 7. Notable features of the three diagrams 
are: (i) the dependence of the maximum value attainable 
by the gas temperature and Mp on the galaxy mass (both 
increase with this); (ii) the existence of a threshold value 
of the metallicity, above which cooling gets so efficient 
that bot temperature and Mp dramatically drop to their 
minimum permitted values, i.e. TcBFt[z(t)] and 0.07M Q , 
respectively; (iii) the dependence of the threshold metallic- 
ity on the galaxy mass (it gets higher at increasing mass). 



Fig. 7. The variation of the peak mass Mp as a function of 
the gas metallicity for all models under consideration. 



6.2. The onset of galactic winds 

In all Type A models the galactic wind occurs at early 
epochs, i.e. within the first Gyr, even though a systematic 
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trend with the galaxy mass can be noticed. High mass 
galaxies undergo their wind ejection earlier than low mass 
galaxies (0.22 Gyr for the 3Ml.t galaxy and 0.82 Gyr for 
the 0.001Ml,t object). In contrast, galactic winds in Type 
B models span a much wider range of ages that go from 
0.22 Gyr for the 3M L . T galaxy to 4.99 Gyr for the 0.1M L>T 
object, to even no ejection at all within the maximum age 
range permitted by the assumed cosmological parameters 
Hq, qo and Zf or as in the case of the two lowest mass 
galaxies. In this latter case, star formation continues all 
over the whole galactic lifetime albeit at a minimal level 
of activity. 

The star formation history of our model galaxies is 
shown in Fig. 8, whereby the solid and dotted lines are 
for Type A and B models, respectively. The rate of star 
formation, in units of M LtRe x 10 12 x M Q /pc 3 /Gyr, is 
plotted against age (in Gyr). For the sake of clarity, we 
limit the age to the first 10 Gyr. With the efficiency of 
star formation we have assumed [v = 10) most of the gas 
is turned into stars within the first few 10 8 years with 
consequent fall-off of the gas content and star formation 
rate. 

The time variation of the gas content is shown in Fig. 9 
where one can notice that even in models with ever lasting 
star formation (the 0.001 and 0.01 Ml,t galaxies of case 
B), the amount of gas left over after the first 0.5 Gyr is 
less than 1% or so. Therefore, in these galaxies the real 
star formation process reduces to very low levels of activ- 
ity. After the minimum in each curve, which corresponds 
to the galactic wind stage, the gas content rises again be- 
cause continuously ejected (into the interstellar medium) 
by dying stars. However, as star formation does no longer 
occur in our models, this gas is supposed to freely escape 
into the intergalactic medium, because it is rapidly heated 
up to high temperatures corresponding to the velocity dis- 
persion of the galaxy (cf. Fabbiano 1986). 

What is the reason for the different behavior between 
models A and B ? and how much do the models of the two 
groups differ each other ? Close inspection of their struc- 
ture shows that there is no significant difference among 
them, the only reason for the difference in t w being the 
wind criterion adopted in Type A models (weaker condi- 
tion) with respect to Type models B (stronger condition). 
It is worth recalling that in Type A models the thermal 
and gravitational energies of the gas are compared to each 
other before cooling, whereas in the Type B ones they are 
compared after cooling. 

It goes without saying that, in Type B models, the 
prolonged activity of star formation will eventually induce 
somewhat high metallicities in the residual gas content. 
This however is less of a problem because only a very 
small fraction of the galaxy mass is stored in stars with 
such high metallicities (see below). 

Although the above uncertainty on t w constitutes a 
point of weakness in our modelling of the interplay be- 
tween gas thermo-dynamics and galactic winds, both the 



results we are going to presents and the very low gas 
content in those galaxies that do not suffer from galac- 
tic winds, hint that this phenomenon is a marginal event 
over the whole evolutionary history of a galaxy not affect- 
ing its global structure. 

It is worth emphasizing that in both cases the relation- 
ship between galaxy mass and age of the galactic wind has 
the opposite trend as compared to the standard supernova 
driven galactic wind model (cf. section 1) in which low 
mass galaxies eject their wind much earlier than the high 
mass ones. We will come back later to this point. 

6.3. What wedge to the energy equation? The role of r\M 

As already anticipated, the term Hm of the energy equa- 
tion (together with the CBR limit) has an important role 
on determining the properties of the models. Why and 
how are easy to understand. They indeed provide a lower 
limit to the temperature attainable by the gas during the 
cooling process. While the effect of TcBii[z(t)] is straight- 
forward, the one of Hm is by far more indirect and dif- 
ficult to understand. This term is constant for a given 
galaxy because so are X and tff are, but it varies from 
galaxy to galaxy. According to its definition, it represents 
a threshold value for the energy input below which gas 
clouds prone to collapse cannot go. Were this term ne- 
glected, owing to the efficiency of cooling (in presence of 
metals, in particular), gas clouds would fall down to ex- 
tremely low temperatures on a very short time scale. Since 
the peak mass of the IMF scales with temperature, Hm 
has a direct effect on Mp and mean slope of the IMF 
above it. The higher Hm, the higher is Mp and the flat- 
ter is the IMF in the mass range above the peak. The 
quantitative effects depend of course on t]m ■ If Hm is effi- 
cient, and Mp in turn is shifted toward the high mass end 
(say above 1M Q ), of the generations of stars born under 
these conditions, a large fraction of these will soon evolve, 
become white dwarf or collapsed remnants in virtue of 
their short lifetime as compared to the evolutionary time 
scale of a galaxy. More gas is refuelled to the interstel- 
lar medium, more metals are created, and star formation 
will tend to last longer. The extreme case can be met in 
which, depending on the competition between gas resti- 
tution and heating, no galactic wind occurs in all model 
galaxies. Therefore Hm eventually drives the mean metal 
content, the relative fraction of living stars and remnants 
(white dwarfs, neutron stars, black holes as appropriate), 
and many other properties of the models, e.g. the colors 
(see below) . All this hints that the parameter t\m can per- 
haps be constrained by looking at global properties either 
of a galaxy, such as the mean and maximum metallicity, 
the colors at the present epoch etc., or of a galaxy man- 
ifold, such as the CMR. Many preliminary models calcu- 
lated to this aim have indicated that t\m = 2 -f- 3 x 10~ 6 
is a plausible choice leading to acceptable models. Other 
values cannot of course be excluded. 
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7. Discussion of the main results 



7.1. Star formation history 

As a result of the complex game between IMF and gas 
thermo-dynamics, the standard supernova driven wind 
model is largely invalidated. Recalling that this latter 
was suggested to explain the CMR as a mass-metallicity 
sequence (the mean metallicity grows with the galactic 
mass), from the entries of Table 1 we see that the same 
trend is now possible even for models in which the overall 
duration of the star forming activity goes inversely propor- 
tional to the galaxy mass i.e. At$F °c Mg 1 as indicated 
by the chemical abundances. Although the age range is 
uncertain because of the above arguments, yet the possi- 
bility exists that, following the dominant initial activity 
of rather short duration, traces of star formation might 
occur over much longer periods of time. Remarkably, this 
is easier to get in low mass galaxies than in high mass 
ones. In this context, an elegant way out can be found 
to the difficulties encountered by standard models as far 
as the distribution of nearby galaxies in the Hp versus 
[MgFe] plane, the trend in the suspected enhancement of 
a-elements, and the gradients in Hp and [MgFe] across 
a galaxy are concerned (see Section 2 and the discussion 
below). We will see that the standard interpretation of the 
CMR is maintained as also these models yield the right 
mass-metallicity sequence. Finally, there is another conse- 
quence of the density dependence of the IMF: since within 
a galaxy the central regions are naturally denser than the 
external ones, we expect over there the same trend as pass- 
ing from a low mass (high density) to a high mass (low 
density) galaxy. Therefore, the central regions of a galaxy 
should form stars for periods of time longer than their 
outskirts. If so, galaxies should be built up by an out-in 
process. 

7.2. Luminous versus dark material 

Since in the course of evolution, the peak mass of the 
IMF first grows to a maximum and then declines, and the 
maximum value can be as high as 1M©, we expect dif- 
ferent proportions of living stars and collapsed remnants 
as a function of time (the present age in particular). In 
brief, of the many stellar generations formed during the 
activity period AtsFR, those with lifetime shorter than 
Tq — AtsFR, where Tq is the current galactic age, are now 
in collapsed remnants (black holes, neutron stars, white 
dwarfs as appropriate to the initial star mass) . They con- 
stitute the "baryonic dark matter" scarcely contributing, 
if nothing at all, to the total light. It goes without saying 
that the remaining stellar content is in form of visible stars 
almost fully contributing to the total light. How much 
mass is present in " baryonic dark matter" ? How does the 
fractionary value of this mass vary with the galaxy mass 
(local density and temperature in turn) ? To this aim, in 
Fig. 10 we plot the fractions of visible stars (solid line) 
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Fig. 8. The normalized rate of star formation SFR as a 
function of the age for Type A (thick solid lines) and B 
(thin dashed lines) models. The SFR is expressed in units of 
Ml, Re x 10 12 x M Q /pc 3 /Gyr. For the sake of clarity the age 
is limited to the first 10 Gyr. 
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Fig. 9. The fractionary mass of gas as a function of the age 
for Type A (thick solid lines) and Type B (thin dashed lines) 
models. For the sake of clarity the age is limited to the first 10 
Gyr. 
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Table 2. Amounts of mass in form of oxygen and iron ejected into the intra-cluster medium by galaxies of different mass Ml,t 
over their whole life. The data refer to Type A models with nu = 2 x 10 -6 . The symbol gw indicates the amount of material 
ejected at the stage of galactic wind. The symbol sic indicates the total amount of material emitted by dying stars from the wind 
epoch to the present time. Fer and Ot are the totals. M| J is the total amount of gas thrown into the intra-cluster medium. All 
masses are in Mq. yFe and yo are crude estimates of the galactic yields of the two elements, i.e. Mj /Mi,,R a . Finally, Ot/Fct 
is a crude estimate of the degree of a-enhacement in the ejected material. 
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Interestingly enough, in the i? e -sphere of a low mass 
galaxy (0.001 Ml,t) about 80% of the original mass Ml, Re 
is in form of visible stars, the rest is in remnants, and there 
are almost no traces of gas. Conversely, in the same region 
of a massive galaxy (3 Ml,t) only 25% of the original mass 
Ml, Re is in now form of living stars, about 20% of it in 
form of collapsed remnants, and 55% of it has remained in 
form of gas partly expelled at the stage of galactic winds 
and partly continuously ejected by dying stars and never 
re-used of form stars. For the same arguments advocated 
to predict the gradient in the duration of the star form- 
ing activity across a galaxy, we expect that the fraction 
of living stars ought to decrease with galactocentric dis- 
tance. It goes without saying that as a consequence of this, 
one expects the mass-luminosity ratio to increase with the 
galaxy mass and the radial distance. 

The time dependence of the gas content deserves some 
final remarks. At the onset of galactic winds, the gas con- 
tent has reached its minimum value because of the past 
star formation activity. However, generous amounts of gas 
are later lost by dying stars. The fractionary mass of the 
newly emitted gas is significantly larger than what cus- 
tomarily possible with standard, Salpeter-like IMF. Since 
this gas is likely soon heated up above the virial tempera- 
ture and lost by the galaxy, it will significantly contribute 
to enrich and fuel the intra-cluster medium (see below. 

7.3. Metallicity and enhancement of a- elements 

How do the metallicity and enhancement of a-elements 
vary as a function of time and galactic mass under this 
continuously changing IMF ? The answer is provided by 
Figs. 11 and 12, in which we plot the relative distribu- 
tion of stars, N[Z(i)], per metallicity bin together with 
the relationships between [Fe/H], [C/Fe], [O/Fe], and 
[Mg/Fe], and the total metallicity Z. Fig. 11 refers to the 




Fig. 10. Fractionary masses in stars (solid lines), and remnants 
(dashed lines) as a function of the age for the 0.001 Ml,t (thin 
lines) and the 3 Ml,t (thick lines) galaxies. 



and collapsed remnants (dashed lines) as a function of the 
age. The corresponding gas fraction was already shown in 
Fig. 9. For the sake of clarity, we limit ourselves to display 
only the 0.001 and 3 Ml,t galaxies. 

It is worth noticing that while the fractionary masses 
of living stars and gas are calculated with high precision 
from the set of equations governing chemical evolution, 
the fractionary mass of remnants is calculated separately 
with a precision of about 5%. This explains why the sum 
of the three fractions may sometimes be slightly greater 
than one (1 -j- 1.05 in general). 
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Table 3. Normalized integrated magnitudes rrii of Type A model galaxies with different ages. Real magnitudes are given by 
Mi — — 2.5 log(Mi i H e x 10 12 ) +m,i. Magnitudes are for the Johnson photometric system. Ml,t and Ml,b, c are the mass of the 
galaxy and the mass inside the effective radius, both are in units of 1O 12 M0; M s and M r are the fractionary masses in living 
stars and collapsed remnants, respectively. The normalization mass is Ml,r c . 
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0.001 Ml,t galaxy, whereas Fig. 12 refers to the 3 Ml.t 
object. Both galaxies belong to group Type A and are 
displayed at the age of 15 Gyr. Similar considerations and 
conclusions apply to Type B models as well. There are 
two important results out of the data shown in Figs. 11 
and 12: 

— In the 3Mij galaxy there are very few stars stored in 
low metallicity bins say up to Z = 0.01. This implies 
that for this galaxy the so-called G-Dwarf analog does 
not occur. In the opposite extreme, i.e. the O.OOlMj,^ 
object, about 50% of the stars are in the bins up to 
Z = 0.01, but [Fe/H] increases rapidly to a tenth of so- 
lar, so that very few stars of very low Z are found. Also 



in this case the G-Dwarf analog is avoided. The kind 
of IMF we are using naturally secures a lower metal- 
licity cut-off in the distribution of living stars, similar 
to what required by BCF94 to remove low metallicity 
stars, and a sort of prompt enrichment in the closed- 
box model. The need of inf all- schemes to avoid the G- 
Dwarf problem is somehow ruled out by the new IMF. 

The totality of stars in the SMl^t galaxy have abun- 
dance ratios [C/Fe] and [O/Fe] well above solar. As 
far as the ratio [Mg/Fe] is concerned, it seems that 
about 50% of stars (those with Z < 0.04) have the ratio 
[Mg/Fe] above solar whereas the remaining stars have 
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Fig. 11. The relative number distribution of living stars 
N[Z(t)] per metallicity bin (solid line). N[Z(t)] is on loga- 
rithmic scale and in arbitrary units. The relationship between 
[Fe/H], [C/Fe], [O/Fe], [Mg/Fe], and the total metallicity 
Z. The galaxy on display is the 0.001 M^t model of Type 
A at the age of 15 Gyr. The largest fraction of stars in this 
galaxy have abundance ratios well below solar, this in spite of 
the short duration of the star forming period (0.82 Gyr). This 
model is the prototype of a galaxy without enhancement of the 
a-elements. 



Fig. 12. The relative number distribution of living stars 
N[Z(t)] per metallicity bin (solid line). N[Z(t)] is on loga- 
rithmic scale and in arbitrary units. The relationship between 
[Fe/H], [C/Fe], [O/Fe], [Mg/Fe], and the total metallicity Z. 
The galaxy on display is the 3 Ml, Re model of Type A at the 
age of 15 Gyr. The largest fraction of stars in this galaxy have 
abundance ratios well above solar. This model is the prototype 
of a galaxy with enhancement of the a-elements. 



[Mg/Fe] below solar. Before drawing any conclusion 
out of this, it should be kept in mind that the stellar 
yield of Mg used by Portinari et al. (1997), TCBMP97, 
and in this paper, underestimates the production of 
Mg (sec Portinari et al. 1997 for all details), so that 
the Mg abundance cannot be safely used as an indi- 
cator of a-enhancement. No difficulty exists with the 
yield of O, whereas there is a marginal discrepancy 
in the yield of C. Therefore, adopting [O/Fe] as en- 
hancement indicator, we can conclude that this model 
can be taken as the prototype of an a- enhanced sys- 
tem. The opposite occurs with the 0.001 Ml : t galaxy 
in which almost the totality of stars have [O/Fe] and 
[Mg/Fe] below solar, whereas a large fraction of these 
(the ones with Z > 0.01) have [C/Fe] above solar. 
This model can be taken as the prototype of a galaxy 
with no enhancement in a-elements. Remarkably, this 
is possible in a system in which the past history of 
stellar activity although shorter than in the high mass 
case, still was confined within 0.82 Gyr, too short a 
time for a significant contamination by Type I super- 
novae with the classical Salpeter law. The new IMF 
makes it possible because at the densities typical of 
the 0.001 Ml : t galaxy, the IMF tends to skew toward 
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Fig. 13. The relation between the current gas fraction 
M g /M L>Re and the current metallicity Z(t)/Z Q (Z Q = 0.020) 
for all Type A models with T)m — 2 x 10 -6 . The minimum value 
in M g /ML,Re corresponds to the galactic wind stage. The in- 
crease in Mg/ML,Re afterwards is due to dying stars of the 
previous generations. 
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M=3 M LT 1/-10 L=15 Gyr 



Z. < Z < 2 z„ 



Fig. 14. The theoretical iV(M) for the 0.001 Ml,t model at 
three different ages: 0.4 Gyr (solid), 1 Gyr (dotted), and 15 
Gyr (dashed). 



Fig. 16. The theoretical N(M) for the 3 M l ,t model at the 
age of 15 Gyr. The top solid line is the total N(M), whereas the 
other lines are the partial contributions by different metallicity 
intervals: Z < Zq (dotted), Zq < Z < 2 x Zq (dashed), and 
Z > 2 x Zq (long-dashed). 



M=0.001 M LT v=lO 



r =5 t=15 Gyr 



Z„ < Z < 2 Z, 



Z > 2 Z„ 



Fig. 15. The theoretical N(M) for the 0.001 M l ,t model at 
the age of 15 Gyr. The top solid line is the total N(M), whereas 
the other lines are the partial contributions by different metal- 
licity intervals: Z < Zq (dotted), Zq < Z < 2 x Zq (dashed), 
and Z > 2 x Zq (long-dashed). 



the low mass end, thus favoring the formation of Type 
I supernovae. 

These results are particularly interesting because two 
of the most intriguing and at some extend demanding con- 
straints imposed by the observational data, i.e. the G- 
Dwarf analog and enhancement in a-elements as a func- 
tion of the galaxy luminosity (mass in turn) are ruled out 
by the same token. 

Finally, in view of the discussion below, we present 
the relation between the current gas fraction and current 
metallicity displayed in Fig. 13 (the metallicity is in so- 
lar units). The reversal in the trend shown by each curve 
beyond the minimum value in the gas fraction reached at 
the wind stage corresponds to the further supply by dy- 
ing stars. Two notable features can be pointed out: (i) the 
metallicity is either ever increasing or has a local minimum 
(effect of the gas restitution by stars of the previous gener- 
ations); (ii) the fractionary total amount of gas ejected by 
dying stars is nearly independent of the galaxy mass, going 
from 0.25 for the 3 M L , T galaxy to 0.18 for the 0.001 M LiT 
object (mean value 0.22). Moving from this diagram and 
the entries of Table 1 , one can estimate the total amount of 
gas and metals made available to the intra-cluster medium 
by each galaxy. This estimate is a lower limit because we 
refer only to the i? e -sphere. The amounts of mass in form 
of oxygen and iron ejected by the galaxies (into the intra- 
cluster medium) are given in Table 2, whereby the mean- 
ing of the symbols is obvious. The data refer only to Type 
A models with t)m — 2 x 10~ 6 for the sake of simplicity. 
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The contributions at the stage of galactic wind and all over 
the subsequent history are shown separately. Remarkably 
the contribution from the post wind phase parallels or 
even overwhelms that from the wind stage. The galactic 
yield (yj — MjjM^^e-, j stands for the generic element), 
of Fe varies from yp e ~ 0.005 to 0.008, depending on 
the galaxy mass, whereas that for oxygen increases from 
y Q ~ 0.0005 for the 0.001 M L , T galaxy to 0.018 for the 3 
Ml : t object. How do these results compare with those by 
Gibson & Matteucci (1997), who presented similar calcu- 
lations? Looking at the two model galaxies in common, 
namely the 1 and 0.001 Ml,t objects, and keeping in 
mind that our models refer to the i? e -sphere (therefore 
the amounts of ejecta are under-estimated) whereas those 
by Gibson & Matteucci (1997) are for the whole galaxy, we 
get the following results (GM and CBPT stand for Gibson 
& Matteucci and the present paper, respectively): 

- 1 M L x- (GM/CBPT) g = 2.2, (GM/CBPT) = 0.3 
and (GM/CBPT)fo = 0.4 compared to their standard 
case, and (GM/CBPT) g = 3.8, (GM/CBPT) = 0.3 
and (GM/CBPT)fo = 0.4 compared to their reduced 
binding energy model. 

- 0.001 M L . T : (GM/CBPT) g = 10, (GM/CBPT) = 10 
and (GM/CBPT)f = 1 compared to their standard 
case, and (GM/CBPT) g = 13, (GM/CBPT) = 0.9 
and (GM/CBPT)fo = 0.1 compared to their reduced 
binding energy model. 

We will come back to this topic in the final discussion. 

7.4- How does this IMF vary from galaxy to galaxy ? 

Owing to the systematic change of the IMF with the phys- 
ical conditions of the interstellar medium in which stars 
are formed, we expect the distribution of stars per mass in- 
terval N(M) to vary with time and hence metallicity and 
from galaxy to galaxy. Figure 14 shows N(M) at three 
different values of the age, namely 0.2, 0.6, and 15.0 Gyr, 
for the O.OOIM^t galaxy. The drop-off along the three 
curves at certain values of the mass corresponds to the 
current value of the turn-off, i.e. 0.98 M Q at 15 Gyr, 2.25 
M at 0.6 Gyr, and 3.55 M & at 0.2 Gyr. Notice how at 
increasing age the part of the N(M) below the turn-off 
gradually steepens. At the age of 15 Gyr the IMF much 
resembles the Salpeter law but for the cut-off below 0.15 
Mq. Figures 15 and 16 shows the inner structure of the 
N(M) as a function of the metallicity for the 0.001M LjT 
and 3Ml,t models, respectively, at the constant age of 15 
Gyr. The solid line is the total N(M), whereas the dot- 
ted, dashed, and long-dashed lines are the partial N(M) 
in the metallicity intervals Z < Zq, Zq < Z < 2 x Zq, and 
Z > 2 x Zq, respectively. Note how the trend is reversed 
passing from the O.OOIMj^t galaxy to the 3Ml,t one. In 
the former, nearly the totality of stars have mctallicitics 
smaller than 2 x Zq , in the latter the opposite occurs with 
the majority of stars in the range Z > Zq. 



8. Photometric synthesis: CMR and FP 

Since the IMF is no longer constant with respect to time 
the standard SSP-technique cannot be used to calculate 
the integrated flux, spectrum, magnitudes, and colors of a 
galaxy (BCF94, TCBF96). They must be derived from di- 
rectly calculating the integrated monochromatic flux gen- 
erated by the stellar content of a galaxy of age T 

F X (T)= / *(t)<t>(t,M)f x (M,r',Z)dtdM (53) 

JO J M L 

where f\(M,r' , Z) is the monochromatic flux of a star 
of mass M, metallicity Z(t), and age t' = T — t, and 
&(t) is the rate of star formation. This poses severe prob- 
lems of numerical nature and computer storage, of which 
no details arc given here. Finally, the spectral library, 
the relationship between effective temperature, metallic- 
ity, and spectral type of the stars in the various evolution- 
ary stages, in usage here are as in BCF94 and TCBF96 to 
whom the reader should refer. 

The photometric data (magnitudes and colors) in sev- 
eral pass-bands of the Johnson system are listed in Table 3 
for three values of the age (15, 8 and 5 Gyr) of Type A 
models. In view of the discussion below, two sets of models 
with r\M = 2 and r\u — 3 are presented. All the magni- 
tudes listed in Table 3 are normalized to the mass M^ne- 
Real absolute magnitudes are therefore given by 

Mi = -2.5log(M L}Re x 10 12 ) + m, (54) 

Since our magnitudes M{ refer to the i? e -sphcre, they 
should be converted to total magnitudes before compar- 
ing them with the observational data. To a first ap- 
proximation, one could apply the maximum correction 
AMj = —0.75 irrespective of the possible of color gra- 
dients. 

8.1. The color-magnitude relation 

The CMR predicted by our models is shown in Fig. 17 for 
the three values of the age and it is compared to the data of 
Bower ct al. (1992) for the Virgo and Coma galaxies. The 
adopted distance modulus to Virgo is (ra — M) = 31.54 
(Branch & Tammann 1992), whereas that to Coma is (to — 
M) = 35.12 (Bower et al. 1992). The thick line labelled 
Data is the linear square fit of the observational data. The 
theoretical CMR are displayed by the linear square fits of 
the data listed in Tables 2 and 3 separately for models 
with T] M = 3 x 10~ 6 (thick lines) and jjm = 2 x 10 -6 
(thin lines). Three values of the age are considered, i.e. 
15, 8, and 5 Gyr. Finally, no correction AMi is applied. 
Although the agreement with the data is not perfect, yet 
the slope of the CMR is matched. Indeed the observational 
data is compatible with the 15 Gyr relationship of the 
rjM = 3 x 10~ 6 models. 

It is worth recalling that even in these models the CMR 
is a mass-metallicity sequence, because despite the shorter 
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Table 4. Integrated magnitudes, colors, and mass to light ratios. 
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duration of the star forming activity and earlier winds in 
massive galaxies with respect to the low mass ones, the 
mean metallicity of the former is higher than in the latter 
thanks to the compensatory effect of the IMF which skews 
toward the high mass end. The great advantage now is 
that the contradiction between CMR and a- enhancement 
as far as the overall duration of the star forming activity 
in galaxies of different mass is concerned (cf. section 2.2 
and 2.3) does no longer occur. 



8. 2. Tilt of the fundamental plane 

The question to be addressed is: can this IMF account for 
the slope of the FP without invoking changes in the virial 
coefficients? To answer this question we have calculated 
the mass-luminosity ratios (M/Lb)q and (M/Lk)q for 
Type A models with r\ M = 2 x 1CT 6 and r\ M = 3 x 10~ 6 . 



and three values of the age (5, 8 and 15 Gyr). All the data 
are listed in Table 4. 

The comparison with the observational data for Virgo 
galaxies is shown in Fig. 18, which displays both the 
log(M I Lb)q versus logMyi plane and the K3 versus K\ 
plane. The filled circles are the data from Bender et 
al. (1992), the thick dashed line is the relation k 3 = 
0.15ki + 0.36 presented in section 2.8, the solid lines are 
the model results for i]m = 3x 10~ 6 , whereas the dotted 
lines are the same but for t]m = 2 x 10 ~ 6 . The FP for 
three values of the age are displayed, i.e. 15, 8, and 5 Gyr. 
Remarkably, there is no appreciable variation of the slope 
of the theoretical relation with the age. 

Assuming that both the theoretical and the observa- 
tional FP are represented by linear relationships of the 
type 



log(M/L B ) = 4x logM 12 + B 



(55) 
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Fig. 17. CMR. The filled and open circles are the data for the 
Virgo and Coma galaxies of Bower et al. (1992). The distance 
modulus to Virgo is (m — M)o = 31.54 (Branch & Tammann 
1992), whereas that to Coma is (m — M)o = 35.12 (Bower et 
al. 1992). The thick solid line labelled Data is the linear square 
fit of the data. The theoretical relationships are displayed by 
the linear square fit of the data of Tables 2 and 3: thick lines 
for models with T)m = 3x 10 -6 and thin lines for models with 
t]m — 2 x 1CP 6 . Three values of the age are considered, namely 
15, 8, and 5 Gyr, as indicated along each line. 

and 

k 3 = a x K\ + /3 (56) 

and casting the theoretical relation in terms of the obser- 
vational one (i.e. a and (3) we get 

log(M/L B )e = ay^ x logM 12 + 7 (57) 
with 

7 - + <*J I x M-) + V3 x /?. (58) 

C2 V 2 c 2 

Therefore, the theoretical slope A corresponds to the ob- 
servational one multiplied by the factor y/3/2, whereas 
the theoretical zero point contains the observational slope 
and zero point, and the virial coefficients C\ and C2 (see 
also Burstein et al. 1997). 

To compare theory with data we calculate the linear- 
square fit of log(M / L b )q and logMxz both for t]m = 
2 x 10~ 6 and ij M = 3 x 10~ 6 separately, and then for 
the two cases lumped together to somehow take the un- 
certainty on this parameter into account. The results 
are summarized in Table 5 which contains the theoretical 



Fig. 18. FP. The theoretical log(M / Lb)q ratios as a function 
of the total mass logM\2 for the ages of 15, 8 and 5 Gyr as 
indicated. The solid lines are for r/M = 3 x 10~ 6 , whereas the 
dotted lines are for T]m — 2 x 10 -6 . The thick dashed line 
labelled Data is the relation for the Virgo galaxies, whereas 
the filled circles are the data by Bender et al. (1992) for the 
same objects. The small dots show a few dwarf ellipticals for 
purposes of comparison. 



slopes and zero points together with their uncertainties, 
and the corresponding slopes in the k 3 — Ki space. The- 
oretical and observational slopes find good agreement for 
the combined sample 2&3xlCP 6 , which means that this 
basic parameter of the models cannot be too different from 
the present values. The scatter in the zero point is about 
±0.04 which is fully compatible with the observational 
a = 0.05. 

What would the results be for r]M either much lower 
or higher than 2-^3 x 10~ 6 ? We recall here that tjm deter- 
mines how low in temperature the cooling process can go 
in a collapsing cloud and fixes Mp, which ultimately drives 
the efficiency of all the other heating sources. Therefore, 
the effects ofrjM are easy to foresee. If r]u << 2-^3x 10 -6 , 
Mp falls down to the lower limit of the IMF, which be- 
comes Salpeter-like, not many remnants are left over and 
the FP flattens out. The effect is more pronounced in 
the high mass (low density) than low mass (high density) 
galaxies. The opposite occurs if r\u » 2-1-3 x 10~ 6 , in 
which case log(M/ Lb)q for massive galaxies gets too high 
and the FP too steep. The many numerical models com- 
puted to check this points show that not only the FP goes 
wrong but also the CMR (because of the effect of ijm via 
Mp on the mean metallicity of a galaxy). 

In addition to r\Mi the CBR limit concurs to set the 
minimum temperature during the collapse of the star 
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forming clouds in case tjm allows for even lower tempera- 
tures. 

On the basis of the above considerations, the tilt of the 
FP and at some extent the CMR as well reflect the thermal 
history of the star forming gas over the early evolutionary 
stages and its variation with the galaxy mass. What really 
matters here is not the parameter i]m per se, but what 
physical mechanism determines the temperature of the gas 
at the star forming stage. We will came back to this below. 



Table 5. Theoretical FP 



A A B a = A x ^/f 



B-Band Age: 15 Gyr 



2 x 1(T 6 0.171 ±0.018 1.026 ±0.061 0.139 

3 x 10~ 6 0.209 ±0.022 1.136 ±0.061 0.170 
2&3xl0~ 6 0.188 ±0.013 1.080 ±0.041 0.153 



K-Band Age: 15 Gyr 



2 x 10~ b 0.004 ±0.008 0.054 ± 0.026 0.003 

3 x 10~ 6 0.007 ±0.002 0.095 ± 0.064 0.005 
2&3xl0~ 6 0.006 ±0.010 0.074 ± 0.031 0.005 



Finally, we have also looked at the log(M j 'Lk)q ratio 
as a function of logM\2'. the slope and zero point for the 
K-passband are given in Table 5. The slope is significantly 
flatter than in the case of the B-passband. 

9. Summary and conclusions 

In this paper we have presented the preliminary investi- 
gation of the effects of Padoan's et al. (1997) IMF, the 
shape of which depends on the physical conditions (tem- 
perature, density, and velocity dispersion) of the medium 
in which stars are formed. The main goal of this study was 
to seek for a picture of galaxy evolution (both chemical 
and spectro-photometric) in which the pattern of proper- 
ties of elliptical galaxies can find a coherent explanation. 
Our major concerns were the CMR, the tilt of the FP, and 
the enhancement of a-elements in relation to the global 
duration of the star formation activity, prior to the onset 
of galactic winds. 

The greatest advantages of this IMF with respect to 
the classical Salpetcr law or similar laws in literature, 
e.g. Scalo (1986 and references therein), Arimoto & Yoshii 



(1987), Kroupa et al. (1993) are: (i) the slope, function of 
the mass range; (ii) the natural mass cut-off below which 
the IMF tends to vanish; (iii) the peak mass Mp that gets 
very small at decreasing temperature and velocity disper- 
sion, and increasing density of the gas (remarkably, un- 
der typical values for these three physical quantities, e.g. 
those holding for molecular clouds in the solar vicinity, 
the Salpeter law is recovered); (iv) the departure of the 
new IMF from the standard ones under different physical 
conditions, such as those likely met in elliptical galaxies 
of different density and velocity dispersion. This brings a 
new leverage to the problems under examination. 

The new IMF has been applied to models of elliptical 
galaxies, in which the radial dependence of mass density 
(both stars and gas) and star formation rate are taken 
into account. The present analysis is limited to the region 
inside R e . However, the results that one would expect in 
other (more external) regions of the galaxies are easy to 
foresee. 

The temperature, density, and velocity dispersion gov- 
erning the IMF have been derived by solving the energy 
equation in which various sources of heating and cooling 
are considered, and by adopting a simple-minded model 
for the thermo-dynamical description of the gas clouds 
prone to collapse and the star formation process. 

Heating is mainly caused by the radiative cooling of 
supernova remnants and stellar winds, whereas that by 
UV radiation from massive stars is found to be marginal 
as nearly all UV flux is re-processed by dust into the far 
infrared (Granato et al. 1997). However, two more sources 
of mechanical nature are considered: the first one [He] is 
active only in the very early stages of galaxy formation 
and evolution and it is meant to take into account that 
part of the energy liberated by the collapse of the primor- 
dial gas that will likely go into heat. It provides a sort of 
initial condition to start with, or in other words the initial 
flame of our models. The second one [Hm] is active at all 
times and somehow takes into account that existing gas 
clouds will mechanically interact with one another thus 
causing additional heating. Both sources of heat would 
require detailed studies that go beyond the scope of this 
paper. Therefore, both heating rates are affect by a certain 
degree of uncertainty and Hm, in particular, contains the 
parameter r\M that has to be chosen in advance and has 
to be determined by comparing theory with observations. 

The temporal evolution of the peak mass over the early 
evolutionary stages drives the whole problem to which the 
minimum temperature attainable by collapsing gas clouds, 
chemical enrichment, cooling, and under suitable circum- 
stances the CBR limit concur. It goes without saying that 
while the quantitative results are expected to depend on 
the details of our models and on the their parameters, in 
particular, the scenario emerging from our analysis does 
not. The following remarks are worth making: 
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(1) History of star formation. The overall duration of 
the star forming activity is inversely proportional to the 
galaxy mass (or directly proportional to the mean density). 
Specifically, At$F oc Mg 1 is shorter in massive galaxies 
and longer in low mass ones. This in spite of the stronger 
gravitational potential in the former, and as a straight con- 
sequence of the varying IMF. The same kind of reasoning 
applies to different regions of a galaxy: star formation ac- 
tivity lasts longer in the center (high density) than in the 
external (low density) regions. 

(2) Galactic winds. The onset of galactic winds occur 
for the same reasons as in the classical scheme, however 
with the opposite trend as far as the galactic mass is con- 
cerned. Early on in massive galaxies, and much later in 
less massive ones. Despite this, the mean and maximum 
metallicities increase with the galaxy mass, thus providing 
the same bottom line for the interpretation of the CMR 
as in the classical scenario (a mass metallicity sequence of 
nearly coeval objects). 

(3) G-Dwarf analog. The excess of low metallicity stars 
is easily avoided independently of the galaxy mass, be- 
cause the relatively top-heavy IMF during the early stages 
secures prompt enrichment of the gas and a scarce popula- 
tion of low-metallicity stars. This makes the infall scheme 
somehow superfluous. 

(4) Enhancement of a-elements. This IMF naturally 
tends to produce different degrees of enhancement at 
varying radial distance, age, and galactic mass. At given 
age and position within a galaxy, the enhancement in a- 
elements increases with galactic mass. This is the sort of 
trend expected from the line strength indices in galaxies 
of different luminosity (mass), cf. section 2.2. At given age 
and galactic mass, the enhancement in a-elements tends to 
be stronger going from the center to more external regions, 
because of the decreasing density favoring more massive 
stars and in turn shorter durations of the star forming pe- 
riods. Apparently this trend opposes to the gradients in 
the line strength indices Mg^ and < Fe > observed in a 
number of elliptical galaxies (cf. Carollo et al. 1993; Car- 
ollo & Danziger 1994a, b). The reality is more complicated 
than this straight conclusion because the line strength 
indices do not simply correlate with the chemical abun- 
dances, but depend also on the age and partition function 
N(Z), number of stars per metallicity bin (cf. Tantalo et 
al. 1997b). It turns out that a region with prolonged star 
formation, higher metal content, and lower chemical en- 
hancement (our central regions) may have Mg^ indices 
stronger than in the more external zones with shorter star 
forming activity, less metal enrichment, and higher chem- 
ical enhancement (cf. Tantalo et al. 1997b for all details). 
Finally, in any case the enhancement in a-elemcnts de- 
creases at increasing age (effect of Type I supernovae) . 

(5) Baryonic dark matter. This IMF tends to favor the 
formation of remnants (black holes, neutron stars, and 



white dwarfs, these latter in particular). The percentage 
of remnants may easily outnumber that of visible stars at 
least at the present epoch. This holds in massive galaxies 
in general and/or in low density regions of all galaxies. 
These stars are obvious candidates to baryonic dark mat- 
ter, with strong physical and cosmological implications. 
Is there any observational hint for a large percentage of 
white dwarfs in the stellar mix of elliptical galaxies ? Bica 
et al. (1996) analyzing the UV excess in a sample of ellipti- 
cal galaxies noticed absorption features centered at about 
1400 A and 1600 A that were common to both strong and 
weak sources. These absorption features characterize the 
spectrum of moderately cool white dwarfs (DA5). Easy ar- 
guments indicate that in order to reach detectability in the 
UV, the relative mass fraction of such white dwarfs with 
respect to the luminous component of the galaxy should 
be of the order of 10:1. 

(6) Tilt of the FP. The tilt and tightness of the FP follow 
from the systematic variation of the IMF with galactic 
mass and the narrow range of temperatures Tp that is 
likely established both by internal and external conditions. 
The observational FP is best matched by old models (say 
15 Gyr) with r\ M — 2 -f- 3 x 10~ 6 and T P in the range 
20 to 40 K. Other values for r\M and Tp in turn would 
predict FPs either too flat or to steep as compared to the 
observational data. 

How a galaxy manages to get 7]m and Tp in the above 
ranges is not clear and perhaps beyond the scope of this 
paper. In the present models Tp results from internal en- 
ergy sources. External sources are not taken into account 
but for the limit set by the CBR. This latter is mostly ef- 
fective in low mass galaxies (IMF skewed toward the low 
mass end) because in the initial stages owing to the scarce 
energy input from supernovae and stellar wind from mas- 
sive stars, the temperature will easily fall below Tcbr- 
For the purposes of the present study, we have adopted 
the Fricdmann model of the Universe with Hubble con- 
stant H = 50 km s -1 Mpc -1 , qo = and red-shift of 
galaxy formation Zf or = 5 (to which an initial value for 
Tcbr — 15 corresponds). Other values for the three pa- 
rameters are obviously possible. They would anyhow lead 
to similar results. The effect of increasing Zf or is straight- 
forward, because at given metallicity (cooling) during the 
early stages, it would be easier for the collapsing clouds 
to fall below TcBR\z{t)}. Needless to say that, increasing 
Zf or , higher values of TcBR.[z{t)] in the very early stages 
are implied. 

Another possibility is that external conditions, for in- 
stance the CBR itself, set the initial gas temperature on 
the top of which internal processes should be added. This 
would lead to the interesting possibility that Tp correlates 
with Zf or . In such a case large variations in Tp over short 
periods of time are to be expected. However, a thorough 
exploration of all possible cases is beyond the aims of this 
paper. 



26 



IMF: new scenario 



Finally, our analysis of the tilt and tightness of the FP 
simply proves that a solution is possible without invoking 
changes in the virial coefficients c\, and c 2 (see section 2.6) 
or other effects of dynamical nature (cf. Ciotti et al. 1996). 
If deviations from virial conditions and/or dynamical ef- 
fects do actually occur, they will be expected to constitute 
a minor correction to the gross trend established by the 
varying IMF. 

(7) On the iron-discrepancy. Although it is beyond 
the scope of this paper to address the question of the iron 
discrepancy raised by the ASCA data, we would like to 
briefly comment that part of the disagreement could orig- 
inate from the assignment of the iron abundance to the 
stellar content of these galaxies. In fact, it stems from 
the adoption of a particular [Fe/H]-Mg2 calibration, in 
which no dependence on age and possible enhancement of 
a-elements is taken into account. We start noticing that 
the general relationship linking the total metallicity Z to 
the iron content [Fe/H] is 



Fe 

f— 1 
[ H 1 



log[^-} - log[^-} + 



-0.8 x 



a 



Fe 1 



0.05 x {-f 



(59) 



in which X Q = 0.707 and Z Q = 0.018 are the adopted 
solar values, and the effect of a-elements enhancement is 
also brought into evidence. This relation reduces to the 
standard one given by Bertelli et al. (1994) when [a/Fe] = 
and X Q = 0.700 and Z & = 0.020 are adopted. 

Considering that in the course of galactic evolution 
both helium and heavy elements should have increased 
according to the popular enrichment law AY/ AZ = 2-^3 
(cf. Pagel et al. 1992), the proper correlation between X 
and Z of the theoretical SSP and the iron content [Fe/H] 
can be established. A recent calibration of the Mg 2 index 
as a function of the age and metallicity Z (and X) for 
SSPs is by BCT96 from whom we derive the following 
analytical approximation 

Fe 

[—} = [7.56 - 0.340 x t + 0.0148 x t 2 ] x Mg 2 + 
H 

-[1.325 + 0.0031 x t + 0.00074 x t 2 ] (60) 

where t is the age in Gyr, and no effect of a-elements is yet 
considered. The above relation holds for 0.008 < Z < 0.05, 
0.20 < Mg 2 < 0.35, and ages older than about 5 Gyr, 
which is fully adequate to our purposes. 

For the age of 15 Gyr, this relation is similar to that of 
Buzzoni et al. (1992) used by Arimoto et al. (1997) in their 
analysis. Since elliptical galaxies are likely to be rich in a- 
elements (cf. section 2.2), before assigning [Fe/H] to the 
stellar content on the basis of its Mg 2 index, the proper 
correction for over-abundance by a-elements should be ap- 
plied, i.e. the value [Fe/H] read off from eq.(60) should 
be decreased by the quantity — 0.8[a/Fe] - 0.05[a/Fe] 2 . 



Considering that [a/Fe] may vary from 0.5 to 1 accord- 
ing to current data (cf. Matteucci 1997), the expected 
decrease in [Fe/H] goes from 3 to 5, which would re- 
duce if not eliminate the disagreement. Looking at the 
iron over-abundance (with respect to the solar value) ver- 
sus log{Lx/ Lb) diagram of Arimoto et al. (1997), this is 
the case for NGC 4406, NGC 4636, and NGC 4472 (the 
galaxies with the highest Lx/Lg ratio). Problems re- 
main with NGC 4374 for which the stellar over-abundance 
of Fe/Fe Q — 0.85 is derived as compared to the value 
Fe/Fe & — 0.11 indicated by ASCA. Only very high over- 
abundances of a-elements would be able to remove the dis- 
crepancy. That the [a/Fe] can vary from galaxy to galaxy 
is indicated by extant data. That this applies to NGC 4374 
(much higher value in particular) is not granted. Passing 
we note that the fine structure parameter of this galaxy 
is 2.3 (Schweizer & Seitzer 1992). The galaxy has suffered 
from a certain degree of dynamical youth or internal reju- 
venation (this is also suggested by its location in Hp ver- 
sus [MgFe] plane, cf. BCT96). How this would affect the 
Mg 2 index is not easily assessable. The problems raised 
by Arimoto et al. (1997) still persist. In any case, we like 
to point out that enhancement of a-elemcnts is a natu- 
ral side product of the new IMF, in particular in the low 
density regions of a galaxy. 

(8) On the intra-cluster gas. The data presented in sec- 
tion 7.3 and Table 2 can perhaps alleviate the difficulties 
encountered by the standard models (cf. Matteucci 1997 
for a recent review of the subject) in simultaneously ac- 
count for amount of iron and gas observed in the ICM. The 
i? e -sphere of the 1 and 3 Ml : t galaxies eject amounts of 
gas, oxigen, and iron going from 1.6 to 6.9xlO n M0, 5.6 to 
2.3 xl0 lo M© and 3.8 to 7.9 xl0 9 M Q , respectively. These 
numbers can be easily doubled considering the contribu- 
tion from the remaining half of the galaxy not taken into 
account here. As compared to the estimates by Gibson 
& Matteucci (1995) for their standard model referred to 
the whole galaxy, for which they get ~ 3.7 x 1O 11 M of 
gas, ~ 1.6 x 10 9 M© of Fe, and ~ 1.7 x 10 10 M Q of oxygen, 
there is a mean increase by a factor of 1.2 to 2.3 in gas, 2.3 
to 4.7 in O, and 3.6 to 7.3 in Fe. Whether this data can 
completely rule out the above difficulty cannot be said 
without detailed model calculations convolved with the 
Schechter (1976) luminosity functions of galaxies. Work is 
in progress to improve upon the present models and cast 
light on this topic (Chiosi et al. 1997). 

(9) Going toward complete models. The models we 
have presented are limited to the region within the effec- 
tive radius, and therefore they are not yet able to properly 
describe complete galaxies nor to make detailed predic- 
tions for the amount and the chemical composition of the 
gas ejected into the ICM. Furthermore, they cannot pre- 
dict properties specific to the very central regions of galax- 
ies, nor spatial gradients in physical quantities. However, 
as the leading parameter is density, the results are some- 



IMF: new scenario 



27 



how independent from the particular region for which they 
have been derived, and may be extended to other regions 
in other galaxies provided they have the same density. 
What would the result be going further out in radius, i.e. 
to regions of lower density ? The answer is obvious and 
supported by preliminary numerical calculations. More ex- 
ternal regions would follow the same trend as passing from 
a low mass (mean high density) galaxy to a high mass (low 
mean density) one. The lower density would make it easier 
to form massive stars and to anticipate the onset of galac- 
tic winds. The material expelled from the outer regions of 
a galaxy should be more enhanced in a-clcmcnts and in 
more generous relative proportions with respect to what is 
left over as stars. If so, the construction of a galaxy can be 
viewed as a sort of out-inward process with star formation 
lasting progressively longer as we move inward. 

(10) Replacing galaxy masses with densities. All the 

discussion carried out in so far has been made using the 
galactic mass as the parameter ranking the various prop- 
erties instead of the mean density despite the fact that 
density was the underlying basic parameter. Indeed, the 
whole discussion could be re-phrased in terms of this lat- 
ter. This makes the results of our analysis even more gen- 
eral in the sense that they be could applied to sub-units 
of the same system (proto-galaxy) with different initial 
mean density of gas. Let us imagine for the sake of argu- 
ment, that a number of such sub-units can exist within 
the total potential well of the pro-galaxy (dark and bary- 
onic matter), and start forming stars more or less at the 
same time in the remote past. A low density cloud will 
have the peak mass of its IMF skewed toward the high 
mass end (Mp can easily be higher than what found in 
our models with the lowest density) and, therefore, will 
soon experience the conditions for local gas ejection and 
interruption of star formation. If gas loss is too high, the 
newly born stellar system may even disgregate dispersing 
its stars into the global potential well. If not, the stellar 
aggregate will survive. In any case a scarce population of 
stars is left over, the vast majority of which will evolve 
on time scales shorter than the Hubble time becoming 
compact remnants. All of the other stars (low mass ob- 
jects) will remain either as a diffuse medium or a low lu- 
minosity aggregate within the potential well. In contrast, 
a high density cloud that preferentially will be located in 
the centre of the potential well, will have the peak mass 
of its IMF more skewed toward the low mass end, and 
therefore will continue to form stars for longer periods of 
time before meeting the conditions for local gas ejection 
and ceasing star formation. Since high density clouds will 
lose much less gas than the low density ones, they will 
survive as formed stellar aggregates eventually merging 
into a single unit at the center of the galaxy. In the above 
picture of galaxy formation, the bulk of star formation 
has taken place in the remote past. How and whether the 
gas left over by the building process in each sub-unit is 



thrown into the intergalactic medium cannot be assessed 
with this simple scheme. Furthermore, the relative number 
of sub-units and how they eventually coalesce into bigger 
aggregates is a typical problem to be addressed with the 
aid of Tree-SPH simulations of galaxy formation and evo- 
lution in a given cosmological context (cf. section 1) and 
goes beyond the aims of the present study. 
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